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This is a general review on the observations and physics of black hole X-ray binaries and micro- 
quasars, with the emphasize on recent developments in the high energy regime. The focus is put 
on understanding the accretion flows and measuring the parameters of black holes in them. It in- 
cludes mainly two parts: (1) Brief review of several recent review article on this subject; (2) Further 
development on several topics, including black hole spin measurements, hot accretion flows, corona 
formation, state transitions and thermal stability of standard think disk. This is thus not a regular 
bottom-up approach, which I feel not necessary at this stage. Major effort is made in making and 
incorporating from many sources useful plots and illustrations, in order to make this article more 
comprehensible to non-expert readers. In the end I attempt to make a unification scheme on the 
accretion-outflow (wind/jet) connections of all types of accreting BHs of all accretion rates and all 
BH mass scales, and finally provide a brief outlook. 
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I will start by defining what I mean by black hole bina- 
ries (BHBs) and microquasars in this article. I decide to 
restrict myself to only a subclass of BHBs, namely, BH 
X-ray binaries (BHXBs) , since these are the only class of 
BHBs known observationally. I will then simply refer mi- 
croquasars as BHXBs for reasons discussed in Section [H] 

Since many excellent, comprehensive and quite up-to- 
date review articles on BHXBs are readily available in lit- 
erature, I feel it is not necessary to write another bottom- 
up and comprehensive review article on the same subjects 
at this stage. I will thus take quite an unusual approach 
in this article. I will first give some concise guides on sev- 
eral representative review articles [TH1|, with some nec- 
essary updates. I will then focus on the further develop- 
ments on several topics I feel deserve more discussions, 
i.e., BH spin measurements (Section jlVp , hot accretion 
flows (Section W\ , corona formation (Section VI I , state 
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transitions (Section VII) and thermal stability of SSD 
(Section VIII). The emphasis is thus put on understand- 
ing the accretion flows and measuring the parameters of 
black holes in them. Some rather general issues on BH as- 
trophysics, such as what astrophysical BHs are and how 
to identify them observationally, are not discussed here 
but can be found from my recent book chapter entitled 
"Astrophysical Black Holes in the Physical Universe" [S] . 
The usual practice of writing a review article is to end 
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FIG. 1. Illustration of a BHXB and microquasar. X-ray emis- 
sion is produced from the central hot accretion disk. A jet is 
normally observed in the radio band. The companion star 
may produce optical emission lines. The disk wind may be 
observed with UV/X-ray absorption lines. 



by listing some outstanding issues and major unsolved 
problems, and then to propose some possible approaches 
to them. I initially did not do this in the first draft. 
The history of astronomy tells us that major progress is 
almost always made by unexpected discoveries and re- 
search results; unpredictability is an essential nature of 
astronomy. This article is not intended to be read by 
funding agencies or proposal reviewers, so I thought I 
did not have to do it. In astronomy, knowing what has 
happened, but looking and doing it differently are far 
more important and effective than following other peo- 
ple's advises. However, the editors of this book suggested 
me to write a brief "outlook" in the end. I thus did it 
nevertheless. 



II. ACRONYMS AND TERMINOLOGY 

In Table H) I list all acronyms used in this article; most 
of these are quite commonly used in this community. 

A BH binary (BHB) is a gravitationally bound binary 
system in which one of the objects is a stellar mass BH 
with mass from several to tens of solar masses (M0); the 
other object, i.e. its companion, can be either a normal 
star, a white dwarf, or a neutron star (NS). In case a 
binary system consists of two BHs, it is is referred to as 
a binary BH system, which is not covered in this article. 
When the companion in a BHB is a normal star, the gas 
from the star may be accreted to the BH and X-rays are 
produced, as a consequence of the heating by converting 
the gravitational potential energy into the kinetic energy 
of the gas, and a BH X-ray binary (BHXB) is referred to 
as such a binary system, as shown in Figure [l] The possi- 
ble existence of BHXBs was first suggested by Zel'dovich 
& Novikov [SlITj. The first BHXB found is Cygnus X-1 
[5| , now a well-studied system among many others found 
subsequently in the Milky Way and nearby galaxies. 

The terminology of microquasar has some twisting 
in it. Historically it was first referred to the BHXB 
1E1740.7— 2942 in the Galactic center region, because a 
double-sided jet was detected from it, mimicking some 



TABLE I. List of Acronyms. 



Acronym 


Definition 


ADAF 


advection dominated accretion flow 


ADIOS 


advection dominated accretion 




inflow/outflow solution 


AGN 


active galactic nuclei 


BH 


black hole 


BHB 


black hole binary 


BHXB 


black hole X-ray binary 


BL3Q 


broad-line-less luminous quasar 


BLR 


broad-line region 


BP 


Blandford-Payne 


BZ 


Blandford-Znajek 


CDAF 


convection dominated accretion flow 


GF 


continuum fitting 


DIM 


disk instability model 


FRED 


fast rise and exponential decay 


FSRQ 


fiat spectrum radio quasar 


GR 


general relativity 


GRB 


gamma-ray burst 


HAF 


hot accretion flow 


HFQPO 


high frequency quasi-periodic oscillation 


HID 


Hardness-Intensity-Diagram 


HIF 


hot inner flow 


ISGO 


inner-most stable circular orbit 


JDAF 


jet dominated accretion flow 


LFQPO 


low frequency quasi-periodic oscillation 


LHAF 


luminous hot accretion flow 


LLAGN 


low-luminosity AGN 


LMG 


Large Magellanic Cloud 


LMXB 


low-mass X-ray binary 


NDAF 


neutrino dominated accretion flow 


NS 


neutron star 


NSXB 


neutron star X-ray binary 


PDS 


power density spectrum 


PL 


power-law 


QPO 


quasi-periodic oscillation 


rms 


root-mean-squares 


RID 


RMS-Intensity-Diagram 


RQQ 


radio-quiet quasar 


SEAF 


Super-Eddington accretion flow 


SLE 


Shapiro, Lightman & Eardley 


SPL 


steep power-law 


SSD 


Shakura-Sunyaev Disk 


SXE 


soft X-ray excess 


TID 


truncated inner disk 


ULX 


ultra-luminous X-ray source 


WD 


white dwarf 


XRB 


X-ray binary 



quasars with similar radio lobes, which have much larger 
scales [S]. Soon after, superluminal jets are observed 
from a BHXB GRS 1915+105 lOJ, which is also re- 
ferred to as a microquasar. Nevertheless both BHXBs 
are quite unusual compared to many others, thus mi- 
croquasars were considered quite unusual. However it 
became clear that microquasars may be quite common 
among BHXBs, since the discovery of a normal BHXB 
GRO J1655— 40 [TT], whose superluminal jets were ob- 
served "12] and showed some correlations with its X-ray 
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emission jT^]. Subsequently several more BHXBs have 
been observed with superluminal jets. At this stage mi- 
croquasars began to be referred to as BHXBs with rel- 
ativistic jets (with bulk motion of about or larger than 
90% of the speed of light), to distinguish them from NS 
X-ray binaries (NSXBs) that only have mildly relativistic 
jets (with bulk motion of about or smaller than 50% of 
the speed of light) [T3]. However the discovery of rela- 
tivistic jets from a NSXB Circinus X-1 made the situa- 
tion complicated: relativistic jets are no longer uniquely 
linked to BHs [15 . 

Now in retrospect, a microquasar can be literally and 
easily understood as the micro version of a quasar; how- 
ever a quasar may or may not be observed with colli- 
mated jets. A quasar has been already understood as a 
special galaxy centered by an actively accreting super- 
massive BH with a mass from millions to billions of Mq , 
and thus its total light output is dominated by the BH's 
accretion process, in a similar way as in BHXBs. The 
production or lack of relativistic jets may have similar or 
even the same underlying physical mechanisms in BHXBs 
and AGNs, though their surrounding environments may 
modify their observed morphologies [16 . It is therefore 
more natural to simply refer microquasars as BHXBs; in 
the rest of this article, microquasar and BHXB are used 
interchangeably. 

We therefore will focus on BHXBs and thus will not 
discuss binary systems producing gamma-rays and some- 
times radio jets, which are most likely high-mass NSXBs 
and in which jets or pulsars' winds interact with the 
wind of its high-mass companion to produce the observed 
gamma-rays [T7K55] . Such systems show very different 
observational characteristics, e.g. the long (1667 days) 
super-orbital modulation with phase offset (about 280 
days) between its X-ray and radio light curves found in 
LS I-h61°303 [24]. 

III. REVIEW OF REVIEWS 

Here I attempt to review the four recent review articles 
[iHl] I consider most useful to readers. Additional infor- 
mation and updates are provided when necessary. Some 
overlaps exists between these review articles, as expected 
and inevitable in bottom-up review articles. To avoid 
repetitions as much as possible in this article, I thus put 
different emphasizes on different articles, with of course 
my personal tastes and perceptions. 

A. The most recent Science Collection 

The recent collection of perspectives and reviews in 
the Science magazine provides excellent introductions to 
and concise summaries of the current state of our un- 
derstanding of BH physics and astrophysics [H [25H27] . 
To the subjects of this article, the most relevant article 
in this collection is the one by Fender and Belloni enti- 
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FIG. 2. A typical Hardness-Intensity-Diagram (HID) of spec- 
tral evolution of a BHXB, following the A->B-^C-^D-^E-^F 
cycle (top). Steady or transient jets are present during the 
A— >-B or B— >-C— >-D stage. No jets are observed, but hot disk 
winds are ubiquitous during the D— >-E stage. (Figure 2 in 
Ref.Q.) 



tied "Stellar-Mass Black Holes and Ultraluminous X-ray 
Sources", which is focused on the observational charac- 
teristics of BHXBs [T] . In particular, the article provides 
an excellent description of the general picture of spec- 
tral evolutions of BHXBs with the Hardness-Intensity- 
Diagram (HID), which is found to be well correlated with 
the observed jets from BHXBs, as shown in Figure [2] 
Actually this cycle is also well tracked by their flux vari- 
ability, represented by the measured root-mean-squares 
(rms) above its average flux, as shown by the RMS- 
Intensity-Diagram (RID) ^1 of the BHXB GX 339-4 in 
Figure [3j which also shows additional horizontal tracks 
at intermediate intensities. Sometimes, a full HID cycle 
does not go into the soft state at all (Figure |4]), perhaps 
due to a failed outburst 

The basic scenario is as follows. During the initial 
stage of an X-ray outburst of a BHXB (A— >-B), which is 
triggered by a sudden increase of accretion rate onto the 
BH, its spectrum is normally hard and steady jets are 
always observed. After reaching its peak luminosity, its 
spectrum begins to soften in a chaotic way and transient 
jets are normally observed (B— >-C— >-D). After this tran- 
sition, the system calms down with a soft spectrum and 
no jets are present (D— >E). Finally the system returns 
to its quiescent state with a hard spectrum accompanied 
with the reappearance of jets (E— >F). Throughout this 
cycle, the presence of hot accretion disk winds appears 
to be anti-correlated with its spectral hardness and jet 
ejection. This empirical pattern appears to be quite uni- 
versal for all BHXBs with very few exceptions, though 
the underlying physics is still not well understood yet. 
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FIG. 3. The Hardness-Intensity-Diagram (HID) and RMS- 
Intensity-Diagram (RID) of the BHXB GX 339-4. The main 
difference from Figure [2] is that there are additional horizon- 
tal tracks at intermediate intensities, a phenomenon known 
as hysteresis of state transitions, which will be discussed in 
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FIG. 4. The Hardness-Intensity-Diagram (HID) of H1743- 
322. The main difference from Figure [S] is that there is a 
complete track at low intensity from the 2008 outburst, which 
did not go into the soft state at all. (Adapted from Figure 7 
in Ref.[29].) 



Nevertheless putting together the above scenario is a very 
significant progress in this field over the last ten years. 
Due to the conciseness of this article, some other impor- 
tant subjects on BHXBs are not fully discussed and many 
original references are also missing. 

Other articles in this collection PSVET] are less relevant 
to the subjects of this article, but are still quite interest- 
ing to read, within the context of BH astrophysics and 
physics. The only other type of astrophysical BHs known 
to exist in the physical universe are supermassive BHs in 
the center of almost each galaxy. Volonteri f27l concisely 
summarized our current understanding on how they are 
formed and grow over the cosmic time; merging of two 
BHs is a key process here. Thorne [25] focused on what 
happens when BHs merge together to produce gravita- 
tional waves, which might be used as a new laboratory 
for studying gravitational physics and a new window for 



exploring the universe. Witten [23] then explained the 
quantum properties of BHs, in particular the basic ideas 
behind Hawking radiation, which might not be impor- 
tant at all for astrophysical BHs. Nevertheless the un- 
derstanding gained through studying the quantum me- 
chanics of BHs actually plays very important roles in de- 
veloping other theories of physics, such as that in heavy 
ion collisions and high-temperature superconductors, as 
vividly described by Witten [2^ . 

B. The most recent ARA&A article on BHXBs, 
with some updates 

The most recent ARA&A article on BHXBs by Remil- 
lard and McClintock entitled "X-Ray Properties of 
Black-Hole Binaries" [2] (referred to as RM06 hereafter) 
provides the most complete, comprehensive and accurate 
review of BHXBs, which actually covers subjects much 
beyond just the X-ray properties of BHXBs. The In- 
troduction of RM06 highlights the initial theoretical and 
observational studies of BHXBs, followed by brief com- 
ments on several main review articles on BHXBs pre- 
ceding this one and the basic properties of BHs within 
the context of general relativity (GR). All main proper- 
ties of BHXBs known at the time are summarized in the 
first table and figure there. In Table |lTj I compile the 
most updated data on all BHXBs currently known, in- 
cluding three BHXBs outside the Galaxy and not in the 
Large Magellanic Cloud (LMC); the currently available 
spin measurements for these BHs are also included for 
completeness. In Table [H] 

/(Af) = PorbifS/27rG = AlBHsm^i/il + qf, (1) 

where Porb is the orbital period, Kq is the semi- 
amplitude of the velocity curve of the companion star, 
Mbh is BH mass, i is the the orbital inclination angle, 
and q = Mc/Mbh is the mass ratio. In Figure [s] 
I show the updated graphical representation of most of 
the BHXBs hsted in Table HU 

The X-ray properties of BHXBs are characterized by 
their X-ray light curves, timing and spectra. Essentially 
all known BHXBs were discovered initially as bright X- 
ray sources, and the majority of them were detected as 
transient X-ray sources with X-ray all-sky monitors. The 
transient properties of some of BHXBs can be interpreted 
by the disk instability model (DIM) (see Ref.[^ and ref- 
erences therein) , which assumes a constant mass transfer 
rate from the mass donor to the accretion disk. However 
the accretion rate from the disk to the compact object, 
i.e. a white dwarf (WD), a neutron star (NS), or a BH, is 
normally lower than the mass supply rate in the disk, so 
mass is accumulated in the disk. When the accumulated 
mass exceeds a certain critical value, a sudden increase 
of accretion rate results in a nova-like outburst. DIM 
is most successful in explaining an outburst with fast 
rise and exponential decay (FRED). However as shown 
in RM06, many observed X-ray light curves of BHXBs 
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TABLE II. Twenty four confirmed BHXBs and their BH masses and spins. Except those marked with references, all other data 
are taken from Remillard & McClintock (2006) |2] (references therein). 



Coordinate 


Common*" 


-^orb 


f{M) 


Mbh 


at 


a: 


Name 


Name/ Prefix 


(hr) 


(Me) 


(Mo) 


(cJ/GM^) 


(cJ/GM^) 


0422+32 


(GRO J) 


5.1 


1.19+0.02 


3.7-5.0 


-1" 




0620-003 


(A) 


7.8 


2.72+0.06 


6.35-6.85 [301 


Qd 


0.12+0.19 [31] 


1009-45 


(GRS) 


6.8 


3.17+0.12 


3.6-4.7 






1118+480 


(XTE J) 


4.1 


6.1+0.3 


6.5-7.2 






1124-684 


Nova Mus 91 


10.4 


3.01+0.15 


6.5-8.2 


-0.04 




1354-64 


(GS) 


61.1 


5.75+0.30 


- 






1543-475 


(4U) 


26.8 


0.25+0.01 


8.4-10.4 




0.75-0.85 |32| 


1550-564 


(XTE J) 


37.0 


6.86+0.71 


8.4-10.8 




0.06-0.54 [33| 


1650-500 


(XTE J) 


7.7 


2.73+0.56 


- 






1655-40 


(GRO J) 


62.9 


2.73+0.09 


6.0-6.6 


0.93'= 


0.65-0.75 [m 


1659-487 


GX 339-4 


42.1 


5.8+0.5 


- 






1705-250 


Nova Oph 77 


12.5 


4.86+0.13 


5.6-8.3 






1743-322 


(H) 






- 




0.2+0.3 1311 


1819.3-2525 


V4641 Sgr 


67.6 


3.13+0.13 


6.8-7.4 






1859+226 


(XTE J) 


9.2 


7.4+1.1:*= 


7.6-12.0 






1915+105 


(GRS) 


804.0 


9.5+3.0 


10.0-18.0 


0.998 


0.98-1.0 [35] 


1956+350 


Cyg X-1 


134.4 0.244+0.005 


13.8-15.8 [36] 


+0.75-'' 


> 0.95 [37] 


2000+251 


(GS) 


8.3 


5.01+0.12 


7.1-7.8 


0.03 




2023+338 


V404 Cyg 


155.3 


6.08+0.06 


10.1-13.4 


-r 




0538-641 


LMC X-3 


40.9 


2.3+0.3 


5.9-9.2 


-0.03 


< 0.3 I3H1 


0540-697 


LMC X-1 


93.8 


0.13+0.05 


9.4-12.5 ^39] 




0.94-0.99 [101 HI] 


0020+593 


IC 10 X-1 


34.9 


7.64+1.26 


> 20 42, 43] 






0055-377 


NGC 300 X-1 


32.3 


2.6+0.3 


> 10 [4| 






0133+305 


M33 X-7 


82.9 


0.46+0.07 


14.2-17.1 gg 




0.84+0.05 m\ 



a: Reported in the first paper on systematic BH spin measurements ^T; b: Reported in the most recent literature; c; 
Postulated to be extreme retrograde Kerr BH, due to the lack of the thermal disk component above 2 keV (GRO J1719-24 
also belongs to this class) [33; d: Postulated to be non-spinning BH, due to the observed low disk temperature; e: BH mass 
of 7 Mq assumed 47 ; / Based on the inner disk radius decrease by a factor of two from the hard state to the soft state 

transition |47j . 



are far more complicated than just FRED and their re- 
current time scales are also not compatible with DIM. 
Perhaps disk truncation (a subject to be discussed ex- 
tensively more later) and mass transfer instability are 
additional ingredients fSOl. 

X-ray emission from a BHXB is variable at all time 
scales, from the free-fall or Keplerian orbital time scale 
of milliseconds near the BH, to the various oscillations 
(some are even related to GR) in the disk with time scales 
from milliseconds to minutes, to the viscous time scales of 
minutes, and to the various instabilities of different time 
scales. Therefore timing studies of BHXBs can probe 
the geometry and dynamics in BHXBs. However, lack 
of coherent signals, such as that observed from pulsars, 
makes it difficult to unambiguously identify the under- 
lying mechanisms from the detected X-ray variabilities. 
Nevertheless power density spectra (PDS) and rms still 
allow us to make progress in understanding the general 
characteristics of a BHXB, in particular when combined 
with its spectral behaviors, as already discussed briefly 
above. 

The most successful understanding of BHXBs so far 
is the description of their thermal X-ray spectral compo- 
nent by the classical Shakura-Sunyaev Disk (SSD) model 



fST] , After applying GR to SSD, the temperature dis- 
tribution in the disk can be obtained |52) . Applying 
this model to BHXBs, one can even measure the spin of 
their BHs [47] , by assuming that the inner accretion disk 
boundary is the inner-most stable circular orbit (ISCO) 
of the BH (a subject to be discussed more in Section [Tv| ). 
However a power-law (PL) component is almost ubiqui- 
tous in the spectra of BHXBs. The interplay between 
these two components results in various spectral states, 
which are also found to be well correlated with their tim- 
ing properties. 

Historically these spectral states have been named in 
many different ways, reflecting mostly how they were 
identified with the observations available at the times. 
In RM06, three states are defined, which show distinc- 
tively different spectral and timing behaviors, as shown 
in Figure [6] based on the RXTE data on the BHXB 
GRO J1655-40. The thermal state has its X-ray spec- 
trum dominated by the thermal disk component and very 
little variability. The hard state has its X-ray spectrum 
dominated by the a PL component and strong variability. 
The steep power-law (SPL) state is almost a combination 
of the above two states, but the PL is steeper. I will keep 
using these definitions throughout this article for consis- 
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GS 2023+338 GS 1354-64 GX 339-4 4U 1543-47 



FIG. 5. Schematic diagram of the dynamically 

confirmed BHXBs, maintained by Dr. Orosz 
(•http: / / mintaka.sdsu.edu / faculty / orosz / web/). The color 
scale for the 17 objects with low mass companions (i.e. 
stars with masses less than about 3Mq) represents the 
temperature of the star. However, the high mass companions 
in Cyg X-1, LMC X-1, LMC X-3, and M33 X-7 are 
considerably hotter, and are thus not well represented in 
color scale. (Adapted from the plot maintained by Dr. Orosz 
(http:/ /mi ntaka.sdsu.edu/faculty/orosz/web/l and Figure 1 
in Ref.^) 
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FIG. 6. Characteristic spectra and PDS of the BHXB 
GRO J 1655-40 in three different states, namely, steep power- 
law (SPL), thermal and hard states. Please refer to Figure 2 
in RM06 for details. (Figure 2 in Ref.P;.) 



TABLE III. Relations between HID, spectral states, jets and 
wind. 

HID Stage State Transition Jets Winds 
A— 5>B hard->SPL steady no 
B— >-C— >-D SPL— >-soft transient transient? 

D— >-E soft decay no steady 

E— >F soft— >hard steady no 



tency. Please refer to Table 2 in RM06 for quantitative 
descriptions of these three states. When a BHXB is in 
quiescence, i.e. not in an outburst, its spectral shape is 
similar to the hard state spectrum. 

Referring back to Figure [2] RM06 found that the HID 
can be reasonably well understood in terms of transi- 
tions betvifeen the above three states; I summarize this 
in Table |III[ After examining the extensive RXTE data 
on spectral evolutions of six BHXBs, the above conclu- 
sion seems to be valid generally. Thus a coherent pic- 
ture seems to emerge: jets are produced only when the 
PL component in the spectrum is strong and winds are 
present only when the PL component disappears; jets and 
winds appear to be mutually exclusive. This suggests 
that transient winds should appear in the SPL state, 
anti-phased with transient jets, if winds quench jets as 
suggested recently [SSI El] ■ 

The robust link between the PL component and jets 
provides an important clue to the jet production mech- 
anism. The PL component is believed to be produced 
in an optically thin and geometrically thick corona, in 



which hot electrons up-scatter the thermal photons from 
the disk. As discussed in RM06, the inner region of the 
disk appears to be truncated quite far away from the 
BH in the hard state; this should be the primary rea- 
son for the low luminosity of this state, as evidence for 
the existence of BHs in BHXBs. However it is still not 
understood how the corona is formed; I will discuss this 
later in Section IVll 

In the final part, RM06 discussed the exciting possibili- 
ties of using BHXBs as probes of strong gravity, including 
confirming these BHXBs contain true BHs, measuring 
the spins of BHs, relating BH spin to the Penrose pro- 
cess and other phenomena, and finally carrying out tests 
of the Kerr Metric. They also discussed in length how 
to measure the spins of BHs and commented the various 



methods of doing so. In Section IV I will re-visit the BH 
spin measurements in details. 

Lacking of a solid surface, a BHXB is not expected 
to produce and actually never observed to have coherent 
pulses. However quasi-periodic oscillations (QPOs) are 
frequently detected in their X-ray light curves; a QPO is 
defined as a 'bump' feature in the PDS, if Q = ly/Av > 2, 
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where v and lS.v are the peak frequency and FWHM of 
the bump, respectively. The observed QPOs are further 
divided into low frequency (LF) (0.1-30 Hz) and high 
frequency (HF) (> 30 Hz) QPOs. Typical LFQPOs can 
be found as the peak around several Hz in the top and 
bottom PDS of Figure [6| In literature, LFQPOs are fur- 
ther divided into several subclasses and their underly- 
ing mechanisms are far from clear at this stage, though 
many models have been proposed and briefly discussed in 
RM06. HFQPOs have been detected from several sources 
at 40-450 Hz. Of particular interests are their stable na- 
ture when detected, which may be linked to either the 
mass or spin or both of the BH in a BHXB. 

In some cases a 3:2 frequency ratio is found. Al- 
though the 3:2 HFQPO pairs could be interpreted in 
some cpicyclic resonance models ^551 156j . there remain 
serious uncertainties as to whether epicyclic resonance 
could overcome the severe damping forces and emit X- 
rays with sufficient amplitude and coherence to produce 
the HFQPOs. A revised model appHes epicyclic reso- 
nances to the magnetic coupling (MC) of a BH's accretion 
disk to interpret the HFQPOs [ST]. This model naturally 
explains the association of the 3:2 HFQPO pairs with the 
steep power-law states and finds that the severe damp- 
ing can be overcome by transferring energy and angular 
momentum from a spinning BH to the inner disk in the 
MC process. 



C. Two comprehensive and long articles on 
modeling accretion flows in BHXBs 

Page limitations to the above review articles did not 
allow in depth discussions on the detailed processes and 
models of accretion flows in BHXBs, which are respon- 
sible for the above described energy spectra, PDS, state 
transitions, jets, and wind. Here I introduce two com- 
prehensive and long articles just on this, by Done P! and 
Done, Gierlihski & Kubota [Ij. 



1. A beginner's guide 

The first is intended to readers who just start research 
in this field [3] . It started with the basic tools in plotting 
spectra and variability, and then described the basic ideas 
and methods used to infer the inner accretion disk radius 
when a BHXB is in a thermal state, in order to measure 
the BH spin 47 . A useful introduction is made on how 
to make various corrections to account for various effects, 
including color correction, special and general relativis- 
tic effects, starting from the original work [47]. Several 
commonly use fitting models, i.e. DISKBB, BHSPEC 
and KERRBB, in the XSPEC package are also briefly 
introduced. An example is given to demonstrate that the 
expected relation Ldisk « 7]^ [47] agrees with the data 
from the BHXB GX 339-4 P], where Ldisk and Ti„ are 
the disk's total luminosity and temperature at the inner 



disk boundary, respectively. 

The hard state is then briefly touched upon, using 
the Advection Dominated Accretion Flow (ADAF) model 
[55] . ADAF naturally explains the hot corona required 
by the observed PL component, especially at very low 
accretion rate. Evaporation at low accretion rate in the 
inner disk region has been proposed as the mechanism 
producing a geometry with a radially truncated disk and 
a hot inner flow; the latter might be the ADAF [60] . 
This also means that when a significant PL component 
is present in the spectrum, BH spin measurement can- 
not be done with the inferred inner disk boundary (see, 
however, counter evidence discussed in Sections |IV| and 



VI ) . This geometry is considered a paradigm that can ac- 



count for many of the observed diverse phenomena from 
spectral evolution to timing properties. 

A brief, yet interesting discussion is given on scal- 
ing up the above models to Active Galactic Nuclei 
(AGNs), which host actively accreting supermassive BHs 
in the centers of galaxies. Applying the insights learnt 
from BHXBs on their spectral evolution, changing disk- 
temperature with accretion rate and BH mass, and disk- 
jet connections, one might be able to understand many 
phenomena beyond the simple AGN unification scheme, 
in which the different observational appearance is all at- 
tributed to an viewing angle difference. 

The continuum emissions of both the disk thermal and 
PL components are modifled by absorptions and added 
by additional spectral features along the line of sight 
(LOS). Absorptions in neutral media produce various 
photoelectric absorption edges, but in ionized media re- 
sult in both absorption edges and lines. In addition to 
recovering the original X-ray emission of a BHXB, mod- 
eling the absorption features is important in learning the 
physical properties of the absorption media, such as col- 
umn density, ionization and velocity along LOS. Winds 
from BHXBs discussed above are always detected this 
way. Several XSPEC fltting models for various kinds 
of absorption edges and lines are also introduced here, 
i.e. TBABS, ZTHABS, TBVARABS, ZVPHABS, 
TBNEW, ABSORI and WARMABS. Figure [Tj shows 
the calculated absorption structures with ABSORI and 
WARMABS; note that ABSORI does not include ab- 
sorption fines, but WARMABS does. 

X-rays interacting with the surrounding medium can 
ionize it and heat the free electrons up by Compton scat- 
tering; hot electrons can also loose energy by interact- 
ing with lower energy photons. A Compton tempera- 
ture of the plasma is reached when the above two pro- 
cesses reaches an equilibrium. The temperature is deter- 
mined by only the spectral shape of the continuum, so 
the heated plasma can escape as winds from the disk if 
the velocity of the ions exceed that of the escape velocity 
at that radius, which is usually very far away from the 
inner disk boundary. As the continuum luminosity ap- 
proaches the Eddington limit, the radiation pressure re- 
duces the escape velocity substantially in the inner disk 
region such that continuum driven winds can be launched 
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FIG. 8. Ionized reflections from a constant density slab, calcu- 
lated with the XSPEC ATABLE model, REFLIONX.MOD, 
which includes the self-consistent line and recombination con- 
tinuum emission, for different values of ^. The inset shows a 
detailed view of the iron line region. For neutral material, 
around 1/3 of the line photons are scattered in the cool up- 
per layers of the disk before escaping, forming a Compton 
down-scattered shoulder to the line. For highly ionized re- 
flection, the upper layers of the disk are heated to the Comp- 
ton temperature so that the spectral features are broadened. 
(Adapted from Figure 22 in Ref.[3].) 



FIG. 7. Calculated absorption structures of a column den- 
sity A'^H = 10^'^ cm"'^ for different ionization parameter 
^ = L/nr^ , where L is the illuminating luminosity, n is the 
density of the medium, and r is the distance from the radi- 
ation source. Top panel: with the ABSORI model that does 
not include absorption lines. Bottom panel: for ^ = 100 with 
both ABSORI and WARMABS, the latter is based on the XS- 
TAR photo-ionization package thus absorption lines are taken 
into account properly. (Adapted from Figure 17 in Ref. [3].) 



almost everywhere in the disk of a BHXB, forming a ra- 
diation driven wind, which is also called thermal wind; 
alternatively winds may also be driven magnetically in 
the inner disk region, but this is much less understood 
yet. In contrast, in an AGN the peak continuum emis- 
sion is in the UV band, which has a much larger opacity 
than Compton scattering in neutral or weakly ionized 
medium during both photoelectric and line absorptions. 
This means the effective Eddington luminosity is reduced 
by large factors and UV line driven wind is easily pro- 
duced at high velocity. This explains the relatively lower 
velocity (hundreds km/s) and highly ionized winds from 
BHXBs, but much higher velocity (thousands km/s to 
0.2 c) and weakly ionized winds from AGNs. 

Material illuminated by X-rays produce both fluores- 
cence lines and reflection features, which depends on 
both the continuum spectral shape and ionization state 
of the material, as shown in Figure [8] calculated with 



the XSPEC ATABLE model REFLIONX.MOD, which 
includes the self-consistent line and recombination con- 
tinuum emission. Note that for highly ionized reflection, 
the Compton heated upper layers of the disk broadens 
the spectral features; these effects are not included in 
the simpler XSPEC model PEXRIV. Replacing the sta- 
tionary slab by a disk around a BH, both the special and 
general relativistic effects will further smear (broaden) 
the spectral features, as shown in Figure [9] for the iron 
line region with different values of ^, inner disk radius ri„, 
and viewing angle «; these are calculated with the RE- 
FLIONX.MOD and then convolved with KDBLUR in 
the XSPEC package. Other similar XSPEC fitting mod- 
els are DISKLINE, LAOR, and KY. The inner disk 
radius can in principle be determined by modeling the 
observed broad iron line features, that in turn can be 
used to measure BH spins, as will be discussed briefly in 
Section |TV] Figure [TO] shows all components of a broad 
band spectrum of a BHXBs, including interstellar ab- 
sorption. 

We thus have two ways discussed so far to determine 
Till of BHXBs by either fitting its thermal continuum or 
the fluorescent iron line feature. However results are not 
always consistent, obtained with different methods, even 
using the same data by different authors. The author 
believes that r^^ increases at low accrete rate, to make 
room for the corona occupying this space. Nevertheless 
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a fairly good summary on the current conflicts and con- 
frontations is made on this issue, which I will discuss 
further in Section HVl 



Energy (keV) 



FIG. 9. Relativistic smearing of the iron line region. The 
sharp peak shown in solid grey lines are the original features. 
Other colored solid lines have different rin marked in units 
of rg = GM/c? for i — 30° (solid lines); similarly those for 
i — 60° are shown in dotted lines. (Adapted from Figure 25 
in Ref iSJ.) 
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FIG. 10. Illustration of all components of the broad band 
energy spectrum of a BHXBs; interstellar absorption is nor- 
mally prominent at low energies. (Adapted from Figure 3b in 
Ref [HI].) 



2. An expert's handbook 

The second one can be considered as a handbook on 
accretion flows on BHXBs, and really lives up to its sub- 
title "Everything you always wanted to know about ac- 
cretion but were afraid to ask" [1] (DGK07 hereafter). 
Besides its much longer length of 66 pages, the main dif- 
ference from the above review articles is that it is fo- 
cused on confrontations between theories and observa- 
tions and intends to depict a coherence picture of the 
accretion physics in BHXBs. In the following I will sum- 
marize briefly the main points and conclusions reached 
in DGK07; those I have reviewed above and will discuss 
more later will be skipped for brevity. 

The underlying physics of DIM for triggering the 
outbursts in BHXBs discussed above is the hydrogen- 
ionization instability, which produces the so-called "S"- 



curve, as shown in Figure 11 irradiation by the inner 
hot disk to keep the outer disk hot is required to pro- 
duce the slow flux decays, e.g. the exponential decays, 
frequently observed in them. The outer disk radius is 
obviously another key parameter, which is determined 
by the tidal instability in a binary system. This mecha- 
nism can explain why BHXBs with high mass compan- 
ions are all persistently bright, since their outer disks are 
always in the upper branch due to the combination of 
their higher average mass transfer rate and inner disk ir- 
radiation. Similarly it also explains the differences and 
similarities between the light curve properties of neutron 
star X-ray binaries (NSXBs) and BHXBs, since a NS has 
a lower mass. It should be noted that the additional sur- 
face emission from the NS may also help to maintain the 
outer disk hot and stay in the upper stable branch [621. 

The SSD prescription assumes that the stress is pro- 
portional to pressure. Because the gas pressure is Pgas oc 
T but radiation pressure is Prad oc T^, so a small tem- 
perature increase causes a large pressure increase when 
^'rad > ^gas, and thus large stress increase, which in turn 
heats the disk even more. The opacity cannot decrease 
effectively to cool it down, so the disk becomes unsta- 
ble when L > 0.06iEdd, where i^Edd is the Eddington 
luminosity of a BHXB. However in most BHXBs their 
disk emissions appear to be stable up to around 0.5LEdd- 
One way out is to assume that the stress is proportional 
to ^Prad-Pgas c>c T^^^, SO the strcss increases slower in 
the radiation pressure dominated regime. Beyond this 
the disk should become unstable, as evidenced by the 
sometimes "heart-beat" bursts of the super-Eddington 
BHXB GRS 1915-1-105; other BHXBs with i^ax > ^Edd 
(e.g. V404 Cyg and V4641) were not observed with 
such instabilities, perhaps due to the combination of in- 
sufficient observational coverage and sensitivity. How- 
ever similar "heart-beat" bursts were also observed re- 
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FIG. 11. "S"-curve due to the Hydrogen-ionization instability. 
At a radius R, the SSD solution gives kT cx rh^^'^, where T is 
the local disk temperature and m is the local mass accretion 
rate. A positive slope means cooling can balance heating, so 
T changes slowly and the disk is stable. In the unstable mid- 
dle branch, thermally emitted photons are absorbed to ionize 
the hydrogen atoms and are thus trapped in the disk, causing 
a large opacity (a). If the mass supply rate min happens in 
the middle range, then the disk experiences a limit-cycle in- 
stability from (1) to (4): (1) T < 10* K, neutral hydrogen has 
low a, rhin > rh, so the disk is built up and T increases slowly; 

(2) r ~ 10'' K, hydro gen atoms begin to be ionized and so 
has high a, T increases rapidly until hydrogen is fully ionized; 

(3) T > 10^ K, ionized hydrogen has modest a, rhin < rh, so 
the disk is eaten out and T decreases slowly; (4) T ~ 10^ K, 
protons and electrons begin to recombine, T decreases rapidly 
until reaching the lower stable branch again. (Adapted from 
Figure 1 in Ref.l4j.) 



cently from IGR J17091-3624, which is Ukely substan- 
tially sub-Eddington unless it is located much beyond 20 
kpc and/or its BH mass is quite small f63". 

Super-Eddington accretion flow (SEAF) can become 
stable again, if the radiation instability is overcome by 
an optically thick ADAF, i.e. the shm disk model, in 
which the trapped photons in the flow is advected in- 
wards, thus balancing the heating generated by viscosity. 
Strong radiation driven winds can be easily produced; 
this can happen in BHXBs and NSXBs (e.g. Z-sources). 
Evidence exists that truncated inner disk (TID) is quite 
common in SEAF and outflow even dominates over inflow 
in SEAF [55 . 

The observed PL component in BHXBs cannot be ex- 
plained by the SSD-like models and thus requires a hot 
accretion flow (HAF). At low accretion luminosity (e.g. 
the quiescent state) , the HAF may be the advection dom- 
inated accretion flow (ADAF) , the convection dominated 
accretion flow (CDAF), or the advection dominated ac- 
cretion inflow/outflow solution (ADIOS). At higher lumi- 
nosity (e.g. the hard state), the original hot and optically 
thin disk solution (i.e. the SLE solution) is unstable, 
because the electron heating efficiency by the Coulomb 



coupling between protons and electrons is too low. The 
luminous HAF (LHAF), however, has the advection as 
a heating source to electrons, so the heating efSciency 
increases and thus electrons can cool the flow more ef- 
fectively. Outflows can be produced in ADAF/ ADIOS; 
collimated jets can also be produced if magnetic fields are 
involved, so an accretion flow may even be jet dominated 
(i.e. JDAF). 

The interplays between the SSD and HAF may be 
responsible for the observed different states in BHXBs 
discussed above. If the PL component is produced by 
Compton up-scattering, then the combination of the op- 
tical depth r and Ch/Cs, where Ch is the heating power 
in electrons and Cg is the cooling power in seed pho- 
tons, can describe the observed variety of spectra. For 
example, the hard state has Ch/Cs 3> 1, but the thermal 
and SPL state have C^/Cs < 1. The location of Hn is 
proposed to be closely related to Ci^/Cs, as illustrated in 
Figure[T2j As discussed above, the steeper PL in the ther- 
mal and SPL states is mostly non-thermal in nature, so 
non-thermal Comptonization is required. Detailed com- 
parisons with data suggest thermal Comptonization also 
cannot be ignored even in the thermal and SPL states. 

In Figure 12 the hot inner flow (HIE) and patchy 
corona are responsible for thermal and non-thermal 
Comptonization, respectively. When a source transits 
from hard to SPL state, nn decreases, HIE is reduced 
but the patchy corona becomes dominant. The tran- 
sition from SPL to soft state is then marked by the 
disappearance of HIE and significant reduction of the 
patchy corona. The above scenario obviously depends 
on two fundamental assumptions: (1) slab-HIF produces 
the non-thermal PL component; (2) slab-HIF is mostly 
located between the TID and the BH. Both assumptions 
are examined exhaustively based on the existing obser- 
vations and their spectral modeling. Slab-HIF is found 
to be consistent with essentially all data. An alternative 
to slab-HIF is that the PL is produced from the jet base 
and beamed away from the disk; however observations 
suggest the PL emission is quite isotropic, thus in con- 
flict with this. TID is found also consistent with data 
when the PL component becomes important, since the 
observed Ldisk ^ deviates significantly from a linear 
relation, suggesting r-m is larger when the PL component 
becomes important. However some studies showed that 
Tin is unchanged in the initial hard state, if the Compton 
up-scattering process is treated properly in Monte-Carlo 
simulations to recover the lost disk photons [65]; this is- 
sue will be further discussed in Section ITVl 

The above developed HIF/TID model based on the 
spectral evolutions of BHXBs can also be applied to ex- 
plain the majority of time variability PDS of BHXBs 
consistently. The TID acts as a low pass filter as it 
cannot response effectively to variations in the HIE, 
so Tin controls the low frequency break of the PDS: 
i^LFB ^ 0.2(r/6)~3/2(TO/10)-i Hz, where r = n„/rg and 
m = M-Q]^/ Mq] this predicts that i^lfb changes from 
0.03 to 0.2 Hz as observed during transitions from the 
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FIG. 12. Left panel: a selection of spectra in different states taken from the 2005 outburst of GRO J1655-40; state names in 
parentheses are defined in RM06. Right panel: proposed accretion flow structures in these different spectra; note the SSD is 
truncated away from ISCO in the bottom two states. (Adapted from Figure 9 in Ref.[2].) 



hard to SPL and soft states if r decreases from 20 to 
6. The HIF/TID model can also explain the observed 
LFQPO variations with luminosity, if the LFQPOs are 
some kinds of characteristic frequencies related to Tin. 
To some extend, the shape of the PDS can also be ex- 
plained by this model. However the observed variability 
rms ~ flux linear correlation and log-normal distribution 
of fluctuations may need additional ingredient, such as 
the proposed propagation fluctuation model. However 
I point out that such correlations and distributions are 
also observed in the light curves of gamma-ray bursts 
and solar flares [Ml El] , which can be produced by the 
generic self-organized criticality mechanism |66| . Similar 
rms ~ flux correlation has also been found for blazars, 
in agreement with the minijets-in-a-jet statistical model 

The spectral and timing properties of weakly magne- 
tized NSXBs are known to have many similarities and 
differences from BHXBs; this is particularly true for 
the atoll sources that have similar ranges of L/L^^^ to 
BHXBs. The essential distinction is that a NS has a solid 
surface, but a BH does not. Observationally the contin- 
uum spectra of atoll sources can be modeled as composed 
of SSD, PL emission, and blackbody emission from the 
NS surface (or the boundary layer between the TID and 
NS surface) ; the former two components are quite similar 
to BHXBs. The spectral evolutions of these NSXBs are 
thus driven similarly by the combinations of r and Cy^/Cs- 
However the blackbody emission is an additional source 
of Cs, so the PL is not as hard as that in BHXBs and 
Tin variations are less effective in changing C\^/Cs', the 
latter means it is more difficult to find evidence of TID 
from modeling only the spectral evolutions in NSXBs. 
On the other hand, the TID/HIF model in NSXBs can 
produce timing behaviors in the same way as in BHXBs 
discussed above, consistent with observations; additional 
timing behaviors, such as coherent X-ray pulsations and 



kilo-Hz QPOs observed in these NSXBs, are caused by 
the rapid spins of the hard surfaces of the NSs. 

IV. FURTHER DEVELOPMENTS ON BH SPIN 
MEASUREMENTS 

A BH predicted in GR can only possess three param- 
eters, namely, mass, spin and electric charge, known as 
the so-called BH no hair theorem. Even if a BH was born 
with net electric charge, its electric charge can be rapidly 
neutralized by attracting the opposite charge around 
it in any astrophysical setting, because the strength of 
electromagnetic interaction is many orders of magnitude 
stronger than that of gravitational interaction. Therefore 
an astrophysical BH may only have two measurable prop- 
erties, namely, mass and spin, making BHs the simplest 
macroscopic objects in the universe. Practically, only 
Newtonian gravity is needed in measuring the BH mass 
in a binary system. However, GR is needed in measuring 
the BH spin. 

The mass and spin of a BH has different astrophysical 
meanings. Its mass can be used to address the question 
of "How much matter (and energy) has plunged into the 
BH?" . However its spin can be used to address the ques- 
tion of "How did the matter (and energy) plunge into the 
BH?" . This is because matter and energy plunged into a 
BH can carry angular momentum, which is a vector with 
respect to the spin axis of the BH. In order to increase 
the total gravitating mass-energy from Mj with zero spin 
to Mf, the added rest-mass must be [69 1 170 ] 

AM = 3Mi[sin-\Mf/3Mi) - sin-^(l/3)], (2) 

and its final spin becomes 
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where J is the BH's angular momentum. Clearly we have 
a* = 1 when M^/Mi = 6^/^; further accretion simply 
maintains this state [70]. Therefore the required addi- 
tional rest-mass to spin a BH from zero to maximum 
spin is AA/ ~ l.SSMj = 0.75Ajff; this is a lower limit to 
the accreted mass [71 . Figure 13 shows function 
of AM. Ignoring Hawking radiation of a macroscopic 
BH, the only way to extract the energy of a BH and thus 
reducing its gravitating mass is by extracting its spin en- 
ergy. Recently, evidence of BH spin energy extraction to 
power relativistic jets has been found, from the observed 
correlation between the maximum radio luminosity and 
its BH spin of a microquasar [72 [73| ; however the aver- 
age jet power is not correlated with BH spin [T^. This 
indicates that jets may be produced by both Blandford- 
Payne (BP) 75j and Blandford-Znajek (BZ) [76^ mech- 
anisms; but the BZ mechanism is more powerful and re- 
sponsible for producing the peak radio luminosity. This 
provides another possible way to estimate a BH's spin 
[73], similar to a recent proposal of using the peak lumi- 
nosity of the disk emission to estimate a BH's spin W7\. 
However there is so far no independent demonstration of 
validness of either of the above two new methods, which 
are thus not discussed further in this section. 



AM/Mi(a.=0) 
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Figure 14 shows the mass reduction, Am = (A/max — 
Mf)/M,„aTT-^max and Mf are the BH's maximum mass 
and final mass, respectively), as a function of extract- 
ing efficiency e of an extreme Kerr BH. For a IQMq BH 
with e = 1 (Am ^ 0.3), the total extracted energy is 
AmMfC^ ~ 10^^ erg and its total gravitating mass is re- 
duced to Afmax/v^ [7T]; this energy could be sufficient 
to power gamma-ray bursts (GRBs). It is thus plausible 
that supercritical accretion onto a newly born BH may 
spin it up and extract its spin energy to power ultra- 
relativistic jets; multiple spin-up and spin-down cycles 
may also happen during one GRB, if the collapsing ma- 
terial is clumpy. 



As shown in Figure [15) the radius of the ISCO of the 
BH, i?iscOi is a monotonic function of the BH spin [75], 
beyond which radius a test particle will plunge into the 
BH under any perturbation; however, in Newtonian grav- 
ity a stable circular orbit can be found at any radius. It is 
thus reasonable to assume that the accretion disk around 
a BH terminates at this radius, i.e., rin = i?isco- There- 
fore, a* can be inferred if one can measure i?isco in units 
of its gravitational radius Tg = GM/(? . Currently three 
methods have been proposed to measure the BH spin in 
BHXBs, and all these methods rely essentially on measur- 
ing i?isco- In case the radiative efficiency (77 = LjMc?) 
can be measured, a* can also be determined this way, as 
shown in Figure |16[ Actually 77 is a very simple function 
of i?isco, i-e, "(] ~ 1/-RiscOj as shown in Figure 17 
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FIG. 13. BH spin a* vs. accreted rest-mass AM, in units the 
final mass Mf(a, = 1) (soUd line, bottom axis), and in units 
of the initial mass Afi(a, = 0) (dashed, top axis). (Adapted 
slightly from Figure 3 Ref.|71|. 
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FIG. 14. Extracted fractional energy Am = (M^ax — 
Mf)/Mniax as a function of efficiency e of extraction of a BH's 
rotational energy. (Figure 4 in Ref.|71|. 



A. The first black hole spin measurement with 
X-ray spectral continuum fitting 

In 1997, I and my colleagues proposed a method of 
measuring the BH spin in BHXBs. It started when we 
tried to measure the mass of the BH in GRO J1655-40 
[111 by measuring from its X-ray continuum spectral 
fitting and assuming the BH is not spinning; we found 
the mass of the BH is around 4M0 [73]. Coincidentally 
just when this paper was going to press, Orosz and Bai- 
lyn |80j announced their accurate measurement of the 
BH mass of GRO J1655-40, which is significantly larger 
than what we found. A co-author of our paper, Rashid 
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FIG. 15. The radius of the innermost stable circular orbit 
(-Risco) of a BH as a function of the spin parameter (a*) of 
the BH, i.e., the dimensionless angular momentum; a negative 
value of a, represents the case that the angular momentum 
of the disk is opposite to that of the BH, i.e., the disk is in 
a retrograde mode. Therefore a* can be measured by deter- 
mining the inner accretion disk radius, if the inner boundary 
of the disk is the ISCO of the BH. 
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FIG. 16. Radiative efficiency as a function of BH spin pa- 
rameter a,. 



Sunyaev urged me to resolve this apparent discrepancy. 1 
immediately realized that a spinning BH of 7Mq in GRO 
J1655-40 would be consistent with the inferred with X- 
ray data, and added a note in proof in the paper suggest- 
ing this possibility [7S] . This turns out to be the first BH 
spin measurement on record. I then invited two of my 
close friends and collaborators, Wei Cui and Wan Chen, 
to join me to apply this method systematically to other 
BHXBs; a new main conclusion in this work was that the 
first microquasar in the Milky Way, GRS 1915+105, also 
contains a spinning BH 47j. 

At the time DISKBB was the only available fitting 
model in the XSPEC package for determining from 
an observed X-ray continuum spectra, if the disk contin- 
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FIG. 17. The solid line shows numerical curve of the radiative 
efficiency as a function of i?isco- The dashed line shows 
r] ~ l/i?isco, which approximates the numerical curve. 

uum is described by the SSD model: L^isk = 47r(Tr?jTJ^, 
where Tin in DISKBB is the disk temperature at rin 
(i.e. the peak disk temperature in SSD) and Ldisk can be 
calculated from the disk flux (after the correction to ab- 
sorption), the distance to the source and disk inclination 
(an issue to be discussed later). To determine the phys- 
ical inner disk radius from the DISKBB parameter ri„, 
several effects must be considered: (1) electron scattering 
in the disk modifies the observed X-ray spectrum; (2) the 
temperature distribution in the disk is not accurately de- 
scribed by the Newtonian gravity as assumed in SSD; (3) 
the observed temperature distribution is different from 
the locally emitted one; (4) the observed flux is different 
from the locally emitted one. The latter three are all due 
to GR effects [55]. For each of the above effects, we intro- 
duced a correction factor, using the best available knowl- 
edge at the time. Since then, several improvements have 
been made to correct for these effects and this contin- 
uum fitting (CF) method is now quite mature in making 
accurate BH spin measurements, given sufficiently high 
quality X-ray continuum spectral measurements and ac- 
curate system parameters of the observed BHXB. 

B. Fiirther developments and applications of the 
continuum fitting method 

This CF method of measuring BH spin has since been 
applied widely to essentially every BHXB with a well 
measured X-ray continuum spectrum showing a promi- 
nent thermal accretion disk component. In particular, 
this method has been improved and incorporated into the 
widely used X-ray spectral fitting package XSPEC, e.g., 
KERRBB 81J, BHSPEC [82 , and KERRBB2 [35l[37]. 
Both the KERRBB and BHSPEC are relativistic mod- 
els, but they have their own drawbacks and advantages 
(See Ref.[55] for a detailed comparison). KERRBB in- 
cludes all the relativistic effects, but it requires to fix the 
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spectral hardening factor. In contrast, BHSPEC could 
calculate the spectral hardening factor on its own; how- 
ever, it does not include the returning radiation effect, 
which turns out to be an important factor in BH spin de- 
termination in BHXBs. KERRBB2 combines both mod- 
els by generating the spectral hardening factor table from 
BHSPEC and using the table as the input for KER- 
RBB. The research group led by Ramesh Narayan of 
Harvard University, Jeffrey McClintock at Smithsonian 
Astronomical Observatory (SAO), and Ronald Remillard 
of Massachusetts Institute of Technology (MIT) [H3] has 
since applied this method and contributed to most of the 
BH spin measurements available in the community, as 
shown in Table [TTl 

The CP method relies on two fundamental assump- 
tions: (1) The measured r-m is uniquely related to i?isco 
of the BH. (2) There is no or negligible X-ray radiation 
from the plunging matter onto the BH beyond i?isco- 
The latter has been studied with numerical simulations 
that include the full physics of the magnetized flow, which 
predict that a small fraction of the disks total luminosity 
emanates from the plunging region [84) . However, in the 
context of BH spin estimation, it has been found that the 
neglected inner light in the CF method only has a mod- 
est effect, i.e., this bias is less than typical observational 
systematic errors ^SJ [H5] . 

The first assumption above requires that the measured 
Tin remains stable as a BHXB changes its spectral state 
and luminosity. However it was noticed that r-m mea- 
sured is usually much smaller, sometimes even smaller 
than i?isco of a extreme Kerr BH in a prograde orbit, 
when the X-ray spectrum contains a significant hard PL 
component, which is believed to be produced by inverse 
Compton scattering of the thermal disk photons in a hot 
corona. We realized that the inferred smaller Tin could 
be due to the lost thermal disk photons in the scatter- 
ing process. We then investigated this problem and con- 
firmed that the inferred Tin can be made consistent with 
Tin inferred from the thermal disk component dominated 
spectrum, if the scattered photons are recovered prop- 
erly by doing detailed radiative transfer in the corona 
[5S]; the same conclusion was also reached by the Har- 
vard/SAO/MIT group independently without knowing 
our much earlier results [87]. Therefore the method of 
BH spin measurement by X-ray continuum fitting can 
also be applied to some SPL state with strong PL com- 
ponent. The stable nature of the measured r-m is proven 
with the textbook case of LMC X-3, when its X-ray lu- 
minosity varied over more than one order of magnitude 
observed in nearly two decades with many different X-ray 



instruments, as shown in Figure 18 



However, when ordered by the observed disk luminos- 
ity = ^disk/iEddi the mcasurcd shows a clear in- 



creasing trend when /d > 0.3, as shown in Figure 19 



The similar trend has also been found in another BHXB 
GRS 1915-^105 [35] and two NSXBs [Ml HI]. It was 
found that r = r^^Jr^ indeed increases physically when 
?D > 0.3, by comparing the evolution of r as a function 
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FIG. 18. The measured r = Vin/r^ and disk luminosity 
~ ^disk/i/Edd (assuming Mbh = WMq) as functions of 
observation time, for the BHXB LMC X-3. (Slightly adapted 
from Figure 1 in Ref.|88,) 



of over a large range for several BHXBs and NSXBs. 
Using the blackbody surface emissions of the NSs in these 
NSXBs is critical in evaluating any possible disk thick- 
ening due to high luminosity that would block at least 
part of the NS surface emission, as well as determining 
the actual NS mass accretion rate, which turns out to be 
much less than the disk mass accretion rate; this suggests 
that the increased radiation pressure is responsible for 
the increase of r and significant outflow when Id > 0.3 
[64] . The same trend is much more pronounced in the 
super-Eddington accreting ultra-luminous X-ray source 



NGC1313 X-2, as shown in Figure [20] together with the 
data from other BHXBs and NSXBs 90J. However the 
exact value of r obtained this way should be taken with 
caution, since the non-negligible energy advection at high 
accretion rate can modify the disk structure in non-trivial 
ways, thus making the SSD prescription inaccurate in 
this case ^T\ . 

Figure [20] also shows that as lu decreases, r again starts 



to increases. However r increases at higher Zd and with a 
different slope for a NSXB than for a BHXB, which can 
be naturally explained as due to the "propeller" effect 
of the interaction between the NS's magnetosphere and 
its accretion disk [92] and the "no-hair" of the BH [90] . 
Figure [21] shows the radiative efficiencies of various sys- 
tems; BHXBs may have either higher or lower efficiencies 
than NSXBs, because a BH has neither solid surface nor 
magnetic field [5"|. 

It is interesting to compare the BH spin results in our 
first paper [37] and the most recent literature for the same 
BHXBs as listed in Table |llj (1) For GRS 1915-1-105, 
A0620-003 and LMC X-3, both results are fully consis- 
tent; (2) For GRO J1655-40, the original result points 
to an extreme Kerr BH {a^, ~ 0.93) [47], somewhat dif- 
ferent from the most recent result of mildly spinning BH 
(a* = 0.65 — 0.75) [32]. However, the original result was 
obtained using the BH mass of 7Mq, about 10% larger 
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than the currently best estimate that was used to ob- 
tain an updated BH spin in the most recent literature 
[32] ■ From Figure 15 it can be seen that a* would be de- 



creased from 0.93 to 0.87, if the BH mass is decreased by 
about 10%; actually a* = 0.85 is also allowed in the new 
estimate (3) For Cygnus X-1, our original conclu- 
sion that a* — 0.75 in the high/soft state and a* — —0.75 
in the low/hard state was based on the assumption that 
Tin decreased by a factor of two when the source made 
a transition from its normal low/hard state to the un- 
usual high/soft state [94]. However the more realistic 
constraint is that rin changed by more than a factor of 
1.8-3.2 during the state transition [M], which implies that 
a* switched from |a*| > 0.85, fully compatible with the 
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FIG. 21. The diagonal line shows a 1/r scaling, calibrated to 
take the value of 0.057 when r = 6. The thick black line is for 
a BH accreting systems. The range of r = 1 — 9 corresponds 
to the ISCO of a BH with different a, ; the radiative efficiency 
ranges between a few to several tens of percents, far exceeding 
the p-p fusion radiative efficiency taking place in the Sun. The 
case for r > 9 corresponds to TID, whose radiative efficiency 
can be extremely low, because energy is lost into the event 
horizon of the BH. The thin solid black horizontal line is for 
the 10% efficiency when matter hits the surface of a neutron 
star where all gravitational energy is released as radiation. 
The thin solid black diagonal line above the point marked for 
"Kerr BH" is for a speculated "naked" compact object [5], 
whose surface radius is extremely small and thus the radiative 
efficiency can be extremely high. (Figure 10.7 in Ref.[5] 



latest result of a* > 0.95 in the high/soft state [37] . 

The above would suggest that a very low tempera- 
ture disk component exists in its hard state spectrum, 
if Cygnus X-1 indeed harbors a Kerr BH and the accre- 
tion disk switches from a retrograde mode to a prograde 
mode when it makes the hard-to-soft state transition. 
For a supermassive BH in the center of a galaxy, it is 
well understood that its mass is mostly gained through 
accretion in its AGN phase [^S] . Therefore random accre- 
tion (between prograde and retrograde modes) tends to 
make its final BH spin close to zero, regardless its initial 
BH spin. For Cygnus X-1, its current BH mass of around 
10-20 Mq cannot be gained through post-formation ac- 
cretion, since its observed average mass accretion rate 
M ~ 2 X 10~^Mq /yr and the age of its companion is 
much less that 10® yr; actually the age of the companion 
is estimated to range from 4.6 and 7.8 million years |96) . 
This means that its post-formation BH mass growth is 
much less than a fraction of its current mass and thus 
its current BH spin must be quite close to that at birth. 
Even if accreting at Eddington rate, to grow its BH spin 
from and the final mass to be the current observed 
value, the timescale is around 3.1 x 10^ yr and the ac- 
creted mass is roughly 7.3 Mq [3^. If its accretions al- 
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ternates between prograde and retrograde modes, then 
its current BH spin should be even closer to its initial 
spin than that in low mass BHXBs that only stay in one 
accretion mode due to roche-lobe overflows. Therefore 
the current high BH spin in Cyg X-1 must be natal; this 
conclusion is also true for the other highly spinning BHs 
in other BHXBs 

In summary, BH spin in BHXBs can now be measured 
reliably with the CF method, when the luminosity of a 
BHXB is between ~ 0.02 and ~ 0.3 in Eddington unit 
and their system parameters are well-known a prior. The 
observed thermal disk spectrum can be modeled directly 
to obtain the BH spin with the available KERRBB2 
model in XSPEC when the X-ray spectrum in dominated 
by this component, i.e., the source is in the thermal state. 
When a significant power-law component is present, the 
inverse Compton scattering process must be taken into 
account to recover the disk photons scattered into this 
PL component, with for example the SIMPL/SIMPLR 
model [371 inZ] in XSPEC. Recently another way to mea- 
sure BH spin using the outburst properties of BHXBs 
has been proposed j77j , which is mentioned briefly in Sec- 
tion |VlTlj however the effectiveness of this method needs 
to be tested. 
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FIG. 22. Phased VIH ellipsoidal light curves corresponding to 
the three states, namely "passive", "loop" and "active" states; 
the data in each state and each wave band are plotted twice for 
clarity. Typical differential photometric error bars are shown 
in the upper right corner of the passive-state light curve for 
each band. To produce purer ellipsoidal light curves, variabil- 
ity on timescales greater than 10 days has been removed from 
loop- and active-state data. (Figure 2 in Ref.[98].') 



C. Uncertainties of the continuum fitting method 

In spite of the tremendous progress made so far on BH 
spin measurements in BHXBs, most of these BH spin 
measurements suffer from considerable uncertainties. Ac- 
tually the major source of these uncertainties comes pri- 
marily from the uncertainties in their BH masses, accre- 
tion disk inclination angles, and distances. The accurate 
BH mass measurement is required because i?isco must 
be in units of r^ — GMjt? in Figure 
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The disk inch- Figure 22 



nation angle and distance are also required because the 
total luminosity of the disk emission L,x\sV is needed, in 
order to estimate the absolute value of a* . 

So far, all BH masses in X-ray binaries (XRBs) have 
been estimated using the Kepler's 3rd law of stellar mo- 
tion, expressed in the so-called the mass function given 
in Equation ([!]). Since the only direct observables are 
Porb and i^Cj both Mq and i have to be determined in- 
directly in order to obtain the BH mass estimate reliably. 
The companion's mass Mq can be determined relatively 
reliably by the observed spectral type of the companion 
star and i can be estimated by modeling the observed 
ellipsoidal modulation of the companion's optical or in- 
frared light curve. The observed ellipsoidal modulation 
is a consequence of exposing different parts of the pear- 
shaped companion star to the observer at different orbital 
phases (see Figs. [T] and [s]); the pear-shape is caused by 
the tidal force of the compact star, which also heats the 
side of the companion star facing it. For details of BH 
mass estimates using this method, please refer to [5]. 

However model dependence and other uncertainties 
(such as accretion disk contamination) cannot be circum- 



vented completely and thus systematic error may exist 
in determining their system parameters. For example, 
three optical states, namely "passive", "loop" and "ac- 
tive" states, have been identified in the normally called 
"quiescent" state of A0620-003 when its X-ray luminos- 
ity is very low; only during the passive state its optical 
light curve modulation is purely ellipsoidal, i.e., accretion 
disk contamination is completely negligible, as shown in 
98j. This means that considerable system- 



atic errors in determining its inclination angle may occur 
unless only the "passive" state data are used. Unfortu- 
nately previous observations of BHXBs used to determine 
their inclinations did not always occur during the passive 
state, thus systematic errors may be common in previ- 
ous resuhs [iS]. Even for GRO J1655~40, of which all 
previous observations were made during its passive state 
[99] . its BH mass measured with different observations, 
or the same data analyzed by different groups are not 
exactly the same, and even not completely consistent be- 



tween them, as shown in Figure 23 which show a scatter 
of about 20-30% to the estimated BH mass, much larger 
than its statistical error of a few percents. 

It should be pointed out that the inclination i in Equa- 
tion ([ij refers to that of the orbital plane of the binary 
system. Because the accreting material initially carries 
the angular momentum inherited from the companion 
star, the formed disk should be co-planar with the or- 
bital plane of the binary system. However, the BH spin 
in a BHXB cannot be changed significantly by accretion 
[71] (and see discussion above on Cygnus X-1). There- 
fore its spin axis may not coincide with the normal direc- 
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tions in the galactic coordinates, are used to infer their 
distances. For example, the distance of GRO J1655-40 
is commonly taken as 3.2 ± 0.2 kpc, based primarily on 
observed absorption lines and somewhat on the dynamics 
of the observed jets. Critical examinations of all avail- 
able data related to its distance, however, are in favor of 
a much closer distance of less than 2 kpc and more likely 
just 1 kpc 108, 109 . Similar conclusion is also reached to 
the distance of A0620-003, revising its distance from the 
commonly accepted 1.2±0.4 kpc to ~ 0.4 kpc, making its 
possibly the closest BHXB known so far |109| : however, a 
distance of 1.06 ± 0.12 pc was preferred in a more recent 
study 30J. Similarly the currently adopted distances of 
many other BHXBs may also have considerable system- 
atic errors. If true, this would change significantly the 
current estimates on their masses and spins. 

Therefore future improvements of the continuum fit- 
ting method depend upon the improved measurements 
on their BH masses, accretion disk inclination, and dis- 
tances. 



D. Future improvements of the continuum fitting 
method 

The mass ratio q can be determined directly according 
to the law of momentum conservation, i.e., 

Mc/A/bh = Kbh/Kc, (4) 



tion of the orbital plane of the binary system. In this case 
the Bardeen-Petterson effect |106| , due to frame-dragging 
by the spinning BH, can rapidly align the normal axis 
of the inner disk region with the spin axis of the BH, 
making the inner disk and binary system mis-aligned. 
Circumstantial evidence already exists for mis-alignment 
between the two axes, because the orbital inclination of 
GRO J1655-40 shown in the upper panel of Figure [23] 
is significantly different from ^ 80° inferred from its rel- 
ativistic jets |107|, if the jets are powered by extract- 
ing the spin energy of the BH via the Blandford-Znajek 
mechanism [75]. Nevertheless, the orbital inclination is 
normally used in place of the disk inclination, which is 
needed in BH spin measurement using the CF method, 
in order to infer the total disk luminosity and calculate 
GR correction factors. However, the inner disk inclina- 
tion is currently not available essentially for any of the 
known BHXBs. For example, a mis-alignment of 10° 
from i — 70° can cause nearly 50% error to the total disk 
luminosity, which will translate into nearly 30% error in 
Tin. A Schwarzschild BH may be estimated to have a* 
falling any where between [—0.5,0.5], if rin is uncertain 



within about 30%, according to Fig.(15) 



Accurate determinations of distances of astrophysical 
objects in the Milky Way are difficult, e.g., for BHXBs 
that are not standard candles. Normally some absorption 
features in their spectra, in conjunction with their posi- 



if the semi-amplitude of the velocity curve of the BH 
Kqh can be observed directly, as illustrated in Figure [l] 
The orbital inclination i can be then calculated using 
Equation ([T]), avoiding any systematics related to the el- 
lipsoidal light curve modeling. 

Since a BH is not directly observable, we can only hope 
to observe any emission or absorption line feature co- 
moving with it. The accretion disk certainly moves with 
the accreting BH. However any line feature of the inner 
accretion disk suffers from the broadening of disk's or- 
bital motion and distortions by relativistic effects around 
the BH, thus making it practically impossible, or diffi- 
cult at least, for detecting the binary orbital motion of 
the BH. Orbital motion of double-peaked disk emission 
lines were observed for NSXB Sco X-1 [110], the BHXBs 
A0620-003 [nmng) and GRS 1124-68 [m]. Unfortu- 
nately a significant phase offset of velocity modulation 
was found from that expected based on the observed 
orbital motion of the companion, though the velocity 
semi-amplitude is consistent with the expected mass ra- 
tio |112j . Soria et al. (1998) observed the orbital motion 
of the double-peaked disk emission line He ii A4686 from 
GRO J1655-40, and found its velocity modulation phase 
and semi-amplitude in agreement with the kinematic and 



dynamical parameters of the system |113j . However one 
major problem in accurately measuring the orbital mo- 
tion of the primary from the observed double-peaked 
emission lines is how to determine reliably the line cen- 
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FIG. 24. Velocity curve of the 39 observed X-ray absorption 
lines with Chandrairom GRO J1655-40, after subtracting the 
line of sight intrinsic velocity of each line at each orbital phase. 
The upper panel marks the velocity of each line with its l-o 
error bar, slightly shifted horizontally for visual clarity. The 
bottom panel shows the weighted average velocity of all lines 
in the upper panel at each phase; the solid curve is the fitted 
velocity curve with its orbital period and phase fixed at the 
values observed previously. (Slightly adapted from Figures 2 
in Ref.pli]) 



ter, because the lines are typically asymmetric and also 
variable. 

We have recently proposed to observe the Doppler 
shifts of the absorption lines of the accretion disk winds 
co-rotating with the BH around its companion star [114j , 
since in many XRBs accretion disk winds are ubiqui- 
tous and appear to be rather stable when observed (e.g., 
in Ref. I115|). We verified this method using Chandra 
and HST high resolution spectroscopic observations of 



GRO J1655-40 (shown in Figure 24 1 and LMC X-3. 



Unfortunately the currently available data only covered 
small portions of their orbital phases and thus do not 
allow better constraints to their system parameters. Fu- 
ture more observations of these two sources and other 
sources with detectable absorption lines from their accre- 
tion disk winds will allow reliable and precise measure- 
ments of the BH masses and orbital inclination angles in 
accreting BHXBs. 

Since the accretion disk has very high density and is 
ionized near the BH, where the majority of the observed 
disk emission is produced due to the very deep gravita- 
tional potential near the BH, scattering of the primary 
disk emission is inevitable. The scattered light is polar- 
ized and its polarization fraction and position angle de- 
pend on the viewing direction (inclination ), scattering 
optical depth and the radius where the scattering occurs 
[116H119] . Ignoring many details, it can be shown that 
the polarization fraction, P(i), of the observed disk pho- 
tons (initial disk emission plus the scattered emission) is 
given by 



1 ^ , cosi 



P(z) 



1 — cos^ i ' 



(5) 



where ^ is a constant depending upon the scattering op- 
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FIG. 25. Polarization fraction of observed accretion disk emis- 
sion as a function of its inclination. 



tical depth; detailed calculations made by Chandrasekhar 
[120] gives P(75°) = 0.04. Note that here the disk 
photons are from the Raleigh-Jeans part of the multi- 
color blackbody spectrum with a characteristic shape of 
f{i') cx v'^/'^, i.e., no GR effect is included. We can there- 
fore find the disk inclination angle by measuring the po- 
larization fraction of this part of the disk emission, as 
shown in Figure [25] At energies above the Raleigh- Jeans 
part of the multi-color blackbody spectrum, the polariza- 
tion is strongly effected by both the inclination and BH 



spin, as shown in Figure 26 The continuum spectra are 
clearly degenerated for the different combinations of in- 
clination and BH spin, but the polarization fraction and 
angle as functions of energy can clearly distinguish be- 
tween them |118j . Therefore X-ray spectra-polarimetry 



observations of BHXBs will certainly make important 
progresses in measuring BH spin. 

Besides using polarization measurements to obtain in- 
ner disk inclination (and BH spin), the broad iron K- 
alpha measurements can also be used to do so |121H123] , 
because iron K-alpha emission is believed to come from 
the fluorescent emission of the disk as a result of illu- 
mination by a hard X-ray source above the disk. Natu- 
rally this emission is sensitive to the disk inclination (and 
BH spin). However, compared to the polarization mea- 
surement method, this method is less straight forward 
and may suffer from systematic uncertainties in deriving 
the disk inclination (and BH spin), because complicated 
modeling of the hard X-ray component and line emissiv- 
ity from the disk is required. Therefore it is essential 
to use all methods discussed above to measure both in- 
clination angles (orbital plane and inner disk) and BH 
spin. Studying the relationship between the results ob- 
tained with different methods is also important in its own 
rights, in order to understand accretion disk physics, its 
interaction with the BH in its center, production of rel- 
ativistic jets, the origin of the BH spin, and ultimately 
the formation mechanisms of BHs and BHXB systems. 

The recent dispute on the distances of GRO J1655-40 
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FIG. 26. Disk continuum flux (a), polarization fraction (b) 
and polarization angle (c) as functions of photon energy for 
different model parameters, i, inclination angle of the inner 
disk in degrees; m, mass accretion rate in units of Eddington 
rate by assuming a 10% conversion efficiency from rest mass 
to radiation in all cases. Other parameters: mass of the BH is 
IOMq; distance to the BH is 10 kpc; spectral hardening factor 
is 1.6. (Adapted from Figures 4-6 and Table 1 in Ref.^llSj) 



and A0620-003 exemplifies the difficulty of determining 
the distances of BHXBs using mostly absorption hues 
|1081 1109j . We have recently suggested a method of us- 
ing the delay time between the X-ray fluxes of an XRB 
and its X-ray scattering halo by interstellar dust to infer 
its distance |124| I125j . However this method may suf- 
fer from our incomplete knowledge of the distribution of 
interstellar medium. Ideally precise astrometry can de- 
termine their distances model-independently, by measur- 
ing their trigonometric parallaxes. Recently the distance 
to Cygnus X~l was determined reliably and accurately 
this way (1.86;l^Q it^ kpc; 126',); which is key to the con- 
sequent measurement of its BH mass and spin [33 ET] . 
Currently it remains challenging to measure the trigono- 
metric parallaxes of objects at distances beyond several 
kpc where most BHXBs are located. Future high pre- 
cision astrometry missions are expected to improve the 
distance estimates to these BHXBs significantly. 

Therefore we expect that the improved measurements 
discussed above on their BH masses, inner disk inclina- 
tion, and distances will allow future improvements of BH 
spin measurement with the CF method. 



E. Possible application of the continuum fitting 
method to AGNs 



More recently, the CF method is also applied to con- 
strain the BH spin in an active galactic nucleus (AGN), 
which is powered by matter accretion onto the central su- 
permassive BH in its center [127] . However the BH spin 
inferred this way for an supermassive BH is quite uncer- 
tain, because: (1) The peak temperature of the accretion 
disk inversely increases with the mass of a BH. Therefore, 
for a supermassive, its temperature would be in the ultra- 
violet energy range which will be strongly absorbed and 
hard to observe; (2) The system parameters of a super- 
massive BH (e.g., the mass of the BH and the inclination 
angle of the accretion disk) have larger uncertainties; (3) 
The uncertain mechanism for some components (e.g., the 
soft X-ray excess) in AGN spectra also increases the diffi- 
culty; (4) Some emission and absorption lines may distort 
the continuum spectrum substantially; and (5) In some 
cases the contribution of its host galaxy to the observed 
total continuum spectrum cannot be removed satisfacto- 
rily. 

The polarized continuum of an AGN should be a pure 
accretion disk continuum, at least in the optical to near 
infrared band J128'. One possible way to measure the 
BH spin in an AGN with the CF method is to combine 
the observed polarized optical to near infrared contin- 
uum spectrum with the observed total UV continuum 
spectrum to get a broad band continuum spectrum of an 
AGN |129j . In principle the broadened Balmer edge fea- 
tures and the total UV spectrum can be used to constrain 
the disk inclination angle and fraction of host galaxy con- 
tamination, respectively [129) . However the quality of 
the currently available data is still insufficient to allow 
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accurate determination of BH spin in AGNs with this 
method. 

Nevertheless the principal method of measuring BH 
spin in AGNs should be using the reflected broad iron 
line and continuum components |121H123"] . The main 
reason is that i can also be determined simultaneously 
and rin obtained this way is already in units of Tg, thus 
avoiding naturally the uncertainties caused by the BH 
mass and inclination in the CF method. However cross 
calibration can be done between the CF and reflection 
fitting methods if both can be applied to the same su- 
permassive BH. 



V. FURTHER DEVELOPMENTS ON HOT 
ACCRETION FLOWS 

In this section I briefly summarize some further de- 
velopments on hot accretion flows, which are believed to 
be responsible for the PL component of the spectra in 
BHXBs. The multi-waveband spectra of the hard state 
of the BHXB XTE J1118-H480 was modeled with the TID 
with hot accretion flow (ADAF) geometry (panel (a) in 
Figure [27| and TID (with ADAF) plus jet model (panel 

— 5x lO^^mn 



(b) in Figure 27), with tod — 0.05 and mjct 
[55] . The steep UV spectrum provides clear evidence for 
a large truncation radius for the SSD (r — 600), and 
the radio to infrared spectrum dominates the jet emis- 
sion, which also contribute to the hard PL component 
[55] . Such low accretion rate gives very low SSD lumi- 
nosity Zd ~ 5 X 10~^ (since the disk radiation efficiency 
r] ~ 1/r; see Figure [Tt] ) , which is far below the turn-over 
luminosity of Id ~ 10~^ shown in Figure 20 Therefore 
the inferred truncation radius of XTE J 1118-1-480 agrees 
with the extrapolation of data points of XTE J1817~330 
down to very low disk luminosity. Large truncation radii 
are also reported from several other sources in the hard 
state (e.g., in Ref. [130]). However, the exact values of 
these truncation radii may have large uncertainties, since 
no direct detection of the inner disk peak emission was 
available, unlike the strong case of XTE Jl 118-^480 [55] . 
For this reason I did not include these reported values in 
Figure [20l 



Panel (b) in Figure 27 also shows that the radiation 
from both the hot accretion flow and the jet contribute 
to the X-ray emission. However the former is roughly 
proportional to m^, whereas the latter to rh. Thus with 
the decrease of to, the contribution from the jet becomes 
more and more important, thus the X-ray radiation will 
be dominated by the jet (when I < 10~^ — 10~^), as 
shown in Figure 28 [55]. The observational data of very 



low-luminosity AGNs clearly show a correlation between 
radio and X-ray with a correlation index of ^ 1.2 |131j . 
in excellent agreement with the prediction shown in Fig- 
ure 
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It is well known that the highest luminosity an ADAF 
can produce is only about 3%LEdd- However the ob- 
served highest luminosity a hard state can reach can 
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FIG. 27. Spectral modeling results for XTE J1118-I-480. (a) 
Outer cool disk plus inner hot accretion flow model, (b) An 
additional jet component is included. rriD = 0.05 and mjet = 
5 X lO'^riiD. (Adapted from Figures 1-2 in Ref. [93] 1 



be 10%LEdd or even higher, which can be described 
by LHAF shown in Figure [29] where several previ- 
ously known solutions of accretion flows are all unified 
in a single scheme [133] . The hard state spectrum of 
XTE J1550-564 with L ~ 6%LEdd is explained by LHAF 
very well, including the X-ray spectral slope and the value 
of the cutoff energy [134] . 



VI. FURTHER DEVELOPMENTS ON CORONA 
FORMATION 

While the formation of the SSD in a BHXB is reason- 
ably well understood, the formation of the corona, which 
is the hot accretion fiow discussed in Section [Vj remains 
less understood. Over the last more than 10 years, Liu 
and her collaborators have developed a model to explain 
the formation and evolution of the corona in a BHXB 
or AGN, in which the mass accretion rate m in units 
of the Eddington ratio drives the variations of the com- 
plex accretion flows by the interaction between the cold 
SSD and the hot corona [60l [T35HI39] . Specifically, the 
coupling between the hot corona and the cold disk leads 
to mass exchange between them. The gas in the thin 
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flow domination - 



FIG. 28. Radio (8.6 GHz)CX-ray (2C11 keV) correlation for 
BHXBs. The observed correlation is shown by the segment 
AB. Segments BCD show the predicted correlation at lower 
luminosities, which approaches that of a pure-jet model, as 
shown by the segment DE. Note that below the point C (~ 
IQ-'^te), the X -ray emission is dominated by the jet and the 
correlation steepens. (Figures 1 in Ref.tl32j) 




FIG. 29. The thermal equilibrium curve of various accretion 
solutions. The accretion rate is in units of MEdd = lOI/Edd/c^ 
and the units of S is g cm~^. The parameters are A/bh = 
IOM0, a = 0.1, and r = 10. (Figures 1 in Ref.[T33]) 



disk is heated up and evaporates into the corona as a 
consequence of thermal conduction from the hot corona, 
or the corona gas condenses into the disk as a result of 
overcooling by, for example, external inverse Compton 
scattering. If m is low, evaporation occurs and can com- 
pletely remove the thin disk, leaving only the hot corona 
in the inner region and a truncated thin disk in outer 
region; this provides a mechanism for ADAF at low to. 
If TO is high, the gas in the corona partially condenses 
to the disk due to strong Compton cooling, resulting in 
disk dominant accretion. The model naturally explains 
the different structures of accretion flow in different spec- 
tral states as shown in Figure [30| [Ml UMi [HTHni] . The 
hysteresis observed in spectral state transitions can also 
be explained by different irradiations from different evo- 
lution history under the same scenario jl45H147] . 
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Figure 30 is significantly different from the illustration 



FIG. 30. A schematic description of the accretion flow 
structures in different spectral states as consequences of disk- 
corona interaction, which is primarily driven by mass accre- 
tion rate dotm. SPL: steep power-law state; VHS: very high 
state; soft: soft state; H/S: high/soft state; hard: hard state; 
L/H: low/hard state; Q: quiescent state. (Adapted from a 
similar figure provided by Prof. Bifang Liu.) 



of accretion flow structures in different spectral state in 
Figure 12 on three aspects: (1) At intermediate to, the 
SSD is broken by ADAF into two parts, an outer disk 
and an inner disk; (2) The inner disk boundary here 
is always located very close to the BH, except in the 
very low to hard or quiescent state, which is very dif- 
ferent from the TID scenario depicted in Figure [l2] (3) 
The corona here covers essentially the whole accretion 
disk, especially the inner disk region, whereas in Fig- 
ure [12] the corona is mostly located inward from the in- 
ner disk boundary. The observed soft X-ray component 
in the low/hard state can be explained by the existence 
of the inner disk [1391 [HOI [1421 [T44] . As I have discussed 
above, the inner disk boundary radius inferred with the 
CF method in the presence of a strong PL component 
is actually consistent with that in the soft state when 
the PL component is weak, after taking into account the 
Compton scattering in the corona [SSI ■ Actually the 
essential assumption behind the broad iron line/reflection 
fitting method of determining BH spin is that the inner 
disk boundary is at the ISCO when the PL component is 
strong. All these tend to support the existence of the in- 
ner disk at intermediate m. However it remains to be 
demonstrated that if the whole SSD is indeed broken 
into the two parts at intermediate to, as illustrated in 
Figure |30j 

Recently it has been suggested that the ADAF and 
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FIG. 31. Solid lines; accretion disk models at different 
accretion rates. Clumpy structures in either ADAF or corona 
(shown as insets) may be developed at intermediate accretion 
rates, corresponding to different spectral states. (Adapted 
from Figures 1 and 2 in Ref. |148| l 



corona shown in Figure [30| may have clumpy structures, 
as shown in Figure |3l] 148 . The "clumpy" model (Fig- 
ure 31 1 has mainly two different consequences from the 
"uniform" model (Figure [30| : (1) The inner disk is tran- 
sient in the "clumpy" model; (2) The "clumpy" model 
can explain the variabilities observed in X-ray binaries 
(such as the state transitions discussed in Section VII) 
and radio-loud AGNs (such as BL Lac objects). 

Determining the structure of the corona in an BHXB 
observationally remains difficult; the fundamental issue 
on whether the corona covers mostly the accretion disk 
or the central compact object still remains unclear so far. 
Since the observed similar spectral and state transitions 
between some NS and BH XRBs are quite similar |149j . 
it is reasonable to assume that they have similar coronae. 
Recently we have used type I X-ray bursts from low-mass 
NSXBs to show that X-ray bursts experience negligible 
Comptonization and that the corona cools rapidly dur- 
ing the rising phase of X-ray bursts and is then heated 
up rapidly during the rising phase of X-ray bursts, as 
shown in the upper panel of Figure 32 for IGR J1747- 
721 |150j . These results suggest that the corona can- 



not cover the central compact object completely (lower 
panel of Figure 32 ) and that the destruction and forma- 



tion time scales of the corona are as short as seconds; 
such short time scales are quite difficult to understand 
in the above discussed evaporation model, in which the 
time scales are related to the viscous time scales of the 
accretion disk. However, this short time scale is consis- 
tent with a corona produced by magnetic reconnections 
in the accretion disk, in a similar way to the solar corona 
heating |151| : this conclusion was based on the inferred 
accretion flow structure of a BHXB shown in panel (A) of 
Figure |33[ in comparison with the atmospheric structure 



of the Sun shown in panel (B) of Figure 33 



It is worth noticing that the purported coronae in both 
XRBs and AGNs have strikingly similar properties, e.g., 
they all have electron temperatures of the order of hun- 




FIG. 32. Upper panel: Anti-correlation between the 
observed type I X-ray bursts from surface of the NS in 
IGR J1747-721 (Adapted from Figures 3 and 2 in Ref. 150 ). 
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FIG. 33. Schematic diagrams of the solar atmosphere (A) 
and accretion disk structure (B). (Figure 2 in Ref. |15l| ). 



dreds keV, in spite that the temperatures, inner disk radii 
and variability time scales of their cold accretion disks 
all scale with their BH masses (and accretion rates) as 
predicted in the SSD model. This means that their coro- 
nae are scale independent. It is perhaps not coinciden- 
tal that the electrons' velocities in a corona are approx- 
imately the same as the Keplerian orbital velocities of 
the inner disk, which are also roughly the same as the 
launching velocities of jets. Of course these velocities are 
also the varialized velocities of the central BHs. It is 
plausible that turbulent small scale magnetic fields lifts 
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FIG. 34. The clumps in the clumpy corona/ADAF of a 
BHXB are channeled into the rotating and wound-up mag- 
netic field lines and then ejected from the system. The 
launched jets should be episodic and have knotted structures. 



the plasma in the "warm layer" shown in Figure [33] to 
form the corona, which is then launched into the jets by 
the rotating large scale magnetic fields through either the 
BP or BZ mechanism; during these processes the mag- 
netic fields are mostly responsible for changing the direc- 
tions of motion of plasmas by the Lorentz force. If the 
corona is clumped as discussed in Section VI [148 , then 



the plasma channeled into the jets should be clumpy and 
thus the ejected jets should be episodic and have knotted 
structures, in agreement with observations; this scenario 
is illustrated in Figure [34] Alternatively episodic ejec- 
tions may also occur in the disk, similar to the coronal 
mass ejection on the Sun [152j . 

Interestingly, the soft X-ray excess (SXE) frequently 
observed in AGNs can be interpreted as the warm layer 
found in XRBs, since the inferred plasmas parameters, 
with electron temperature of 0.1 to 0.3 keV and Comp- 
ton scattering optical depth of around 10, are similar to 
that of XRBs and universal among AGNs with different 
BH masses, but seem to be related only with rh |153H155] . 
Therefore their warm layers are also scale independent. 
Indeed, magnetic flares, perhaps caused by magnetic re- 
connections in the cold disk, can produce the warm layers 
in both XRBs and AGNs [TMlfT?/] . 



VII. FURTHER DEVELOPMENTS ON STATE 
TRANSITIONS 

To further understand the mechanism of the spectral 
state transitions, it is important to study those spectral 
state transitions in individual BHXBs during different 
outbursts. The advantages are obvious - uncertainties in 
our estimates of black hole mass, the binary properties 
such as the orbital period, or the source distance will not 
play any role in producing the observed diverse transi- 



tion properties in individual BHXBs, which can only be 
driven by the accretion process under different initial con- 
ditions. The sources which firstly allowed such a study 
were the BH transients GX 339-4 and XTE J1550-564, 
the NS transient Aquila X-1, and the flaring NS low-mass 
XRB (LMXB) 4U 1705-44, in which a remarkable corre- 
lation between the luminosity of the hard-to-soft transi- 
tion and the peak luminosity of the following soft state 
was found |158H160] . More recently, a comprehensive 
study of the hard-to-soft spectral state transitions de- 
tected in all the bright XRBs in a period of about five 
years with simultaneous X-ray monitoring observations 
with the RXTE/ASM and the Swift /BAT confirmed the 
correlation between the hard-to-soft transition and the 
peak luminosity of the following soft state |161[ 1162) , as 
shown in the upper panel of Figure 35 [ 161| . More im- 
portant was the discovery of the correlation between the 
transition luminosity and the rate-of-change of the lumi- 
nosity during the rising phase of an outburst of transients 
or a flare of persistent sources, as shown in the lower 
panel of Figure 35 |161j . 

The above correlation implies that in most cases it 
is the rate-of-increase of the mass accretion rate, rather 
than the mass accretion rate itself, determines the hard- 
to-soft spectral transition; this is depicted in Figure [36] 
|161j . In addition, the discovery of the relation between 
the hard X-ray peak flux and the waiting time of tran- 
sient outburst in the BH transient GX 339-4 (as shown 



in Figure 37 1 supports that the total mass in the disk de- 



termines the peak soft state luminosity of the following 
outburst |159j . Therefore spectral states should be un- 
derstood in the non-stationary accretion regime, which 
would be described by both rh and m, as well as the 
total mass in the disk before an outburst. 



VIII. FURTHER DEVELOPMENTS ON 
THERMAL STABILITY OF SSD 

The SSD model predicts that when the accretion rate 
is over a small fraction of the Eddington rate, which cor- 
responds to L > 0.06LEdd, the inner region of the disk 
is radiation-pressure-dominated and then both secularly 
[163j thermally unstable jl641 1165j . However, observa- 
tions of the high/soft state of black hole X-ray binaries 
with luminosity well within this regime (O.OlLEdd ^ L < 

0. 5Lsdd) indicate that the disk has very little variability, 

1. e., quite stable |166j . It has been well established that 
the accretion flow in this state is described by the SSD 
model [31 1167j . Radiation magnet ohydro dynamic simula- 
tions of a vertically stratified shearing box have confirmed 
the absence of the thermal instability |168j . Recently, 
the thermal stability is revisited by linear analysis, by 
taking into account the role of magnetic field in the ac- 



cretion flow |169| . By assuming that the field responds 
negatively to a positive temperature perturbation, it was 
found that the threshold of accretion rate above which 
the disk becomes thermally unstable increases signifi- 
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FIG. 36. Correlation between the peak fluxes of the ini- 
tial low/hard states in the outbursts of GX 339-4 and the 
time since the latest low/hard state peak in the previous out- 
burst. The dashed line passes the origin and the data point of 
maximal peak flux, showing an example of a linear relation. 
(Adapted from Figure 3 in Ref.fTH^l 
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FIG. 35. Upper panel: correlation between the transition 
luminosity (15-50 keV) and the peak luminosity of the follow- 
ing soft state (2-12 keV) in Eddington units. Lower panel: 
correlation between the peak luminosity of the soft state and 
the maximum rate-of-increase of the X-ray luminosity around 
the hard-to-soft transition. (Adapted from Figures 24 and 27 
in Ref.fTeTI) 



cantly, compared with the case of not considering the 
role of magnetic field. This accounts for the stability of 
the observed sources with high luminosities, as shown in 
Figure 38 If the magnetic pressure is less than about 



24% of the total pressure, then this model can explain 
the "heart-beat" limit-cycle instability observed in GRS 
1915-1-105 at its highest luminosity; this peculiar source 
holds the highest accretion rate (or luminosity) among 
BHXBs. 

Observations of GRS 1915-1-105 showed that thigh (the 
duration of the outburst phase) is comparable to tiow (du- 
ration of the quiescent phase) and Lhigh is 3 to 20 times 
larger than Liow [170] . However, numerical calculations 



FIG. 37. A schematic picture of the regimes of the hard state. 
Two assumed transient outbursts of different peak luminosi- 
ties are shown. When a source is under stationary accretion, 
spectral transitions between the hard state and the soft state 
occurs at a nearly a constant luminosity Lq. When a source 
is undergoing an outburst or flare, the hard-to-soft transition 
occurs at a luminosity above Lq. The additional luminosity 
roughly proportional to . The soft-to-hard transitions are 
expected to occur around Lq. (Figures 28 in Ref. [16lj ) 



showed that thigh is less than 5 percent of tiow, and Lhigh 
is around two orders of magnitude larger than Liow • Some 
efforts have been made to improve the theory in order 
to explain observations either by some artificial viscos- 
ity prescription |171j or by additional assumption of the 
energy exchange between the disk and corona [172[ . Tak- 
ing into account the stress evolution process, it was found 
that the growth rate of thermally unstable modes can de- 
crease significantly owing to the stress delay, which may 
help to understand the "heart-beat" limit-cycle variabil- 
ity of GRS 1915+105 [T73]. The Hmit-cycle properties are 
found to be dominated by the mass-supply rate (accre- 
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FIG. 38. The thermal equihbrium curves of a thin disk 
at lOrg for different $i, which is defined as "l>i = B^H and 
is in unit of Gs-cm. Other model parameters are BIf mass 
Mbh = IOM0 and viscosity parameter a = Q.l . (Figure 2 in 
Ref.[T69l) 



tion rate at the outer boundary) and the value of the a- 
viscosity parameter in the SSD model that assumes that 
the viscous torque is proportional to the total pressure 
[77]. It was also found that only the maximal outburst 
luminosity (in Eddington units) is positively correlated 
with the spin of a BH, providing another way to probe 
BH spin |77j ; this is mainly due to the smaller inner disk 
radius and thus higher radiative efficiency for larger a* 
as shown in Figures [15] [16] and [17] 



IX. UNIFICATION AND OUTLOOK 

Despite of the many progresses made over the last 
decades on the study of BHXBs and microquasars, there 
are still many outstanding and unresolved issues, which 
can be pierced by two big pictures. One big picture is 
related to the astrophysics of BHXBs and microquasars, 
which is centered on understanding the accretion-outflow 
(wind or jet) connections at different accretion rates for 
different types of astrophysical systems involving BHs. 
The other big picture is related to the fundamental 
physics of BHXBs and microquasars, which is centered 
on identifying astrophysical BHs and testing theories of 
strong gravity. These two pictures are actually also en- 
tangled together, because a strong gravity theory, such 
as GR, is needed to describe the astrophysical aspects 
of these systems, and we need to understand the astro- 
physics of BHXBs and microquasars before we can start 
to test any strong gravity theory [174]. 

Figure [39] shows the relation between the collimation- 
corrected gamma-ray luminosity and the kinetic power 
for AGNs and gamma-ray bursts (GRBs) [175j . For at 
least some of GRBs, a super-Eddington accreting com- 
pact stellar object, probably a BH with around lOM©, is 
believed to be responsible for powering the highly rela- 
tivistic jets, which produces intense gamma-rays through 
violent collisions of blobs in the iets |176j . The two types 
of AGNs shown here, flat-spectrum radio quasars (FS- 
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FIG. 39. The relation between the coUimation-corrected 
gamma-ray luminosity and the kinetic power for AGNs and 
GRBs. The shaded regions display the 2a confidence band 
of the fits. The blazar and GRB best-fit models (dashed 
and dotted lines, respectively) follow correlations which are 
consistent, within the uncertainties, with the best-fit model 
obtained from the joint data set (solid line). (Figure 3 in 
Ref.[T75]) 

RQs) and BL Lacs, have all been observed to have mildly 
relativistic jets. Their central engines are supermassive 
BHs with around 1Q'^~^Mq with accretion rates just be- 
low (for FSRQs) or far below (for BL Lacs) Eddington 
rate |177j . The good correlation over about 10 orders of 
magnitudes between these very different systems with 
very different accretion rates suggests that there must be 
some common mechanisms responsible for the accretion- 
outflow (wind/jet) connections for all accreting BHs |178j . 

Putting together all related phenomenologies and some 
theoretical modeling of BHXBs and microquasars dis- 
cussed in this article, an unification scheme is illustrated 
in Figure [40] with the major elements in the accretion- 
outflow connections in different types of astrophysical 
systems harboring both stellar mass BHs and super- 
massive BHs with accretion rates over several orders of 
magnitudes. The types of accreting stellar mass BHs 
include BHXBs and microquasars in different spectral 
states, as well as ultra-luminous X-ray sources (ULXs) 
(some of which are most likely ultra-Eddington accreting 
BHs jlSOj ). and GRBs. The types of accreting supermas- 
sive mass BHs include low-luminosity AGNs (LLAGNs), 
BL Lac objects, normal radio quiet quasars (RQQs), FS- 
RQs and broad-line-less luminous quasars (BL3Qs) [179] 
(I create the acronym "BL3Q" here just for fun). 

This scheme is not a theoretical model, or even toy 
model yet. It, however, can be used as one possible 
chain to pierce many observed phenomenologies together, 
providing a possible frame work for further theoretical 
developments on accretion-outflow connections. Future 
observations will scrutinize this unification scheme and 
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Some necessary elaborations and points: 



• m ^ 10: the two-phase clumpy thick disk is in the form 
of neutrino dominated accretion flow (NDAF) , which can 
extend to ISCO due to the extremely low opacity of neu- 
trinos. The large scale magnetic field lines are wound-up 
by the spinning BH and rotating disk. The discontinuous 
jet is produced by the BZ mechanism. Most likely no such 
ultra-Eddington AGN exists. 

• m ~ 1: the thin disk is truncated by radiation pressure, 
which drives near-spherical winds out. BL3Q is predicted 
as the early phase of a quasar, when the radiation pressure 
drives all gas in the radiation cone out so that no broad- 
line region (BLR) can be formed fl 79 j . Actually all gas 
mixed with dust in the radiation cone can be blown out 
at substantially below Eddington rate, since the dust has 
a very high opacity to emissions from visible to UV. There 
BL3Q should be generic in the early phase of a quasar's 
active cycle, if the quasar activity is triggered by enhanced 
mass infiow consisting of a mixture of gas and dust. 

• m ~ 0.3: the thin disk extends to around ISCO and the 
two-phase corona is clumpy, so the jet is produced mostly 
via the BZ mechanism and is discontinuous. The dusty 
torus and BLR for in a FSRQ are not shown. 

• m ~ 0.1: the thin disk extends to around ISCO and near- 
equatorial wind is thermally driven out; small scale and 
turbulent magnetic fields may be responsible for launching 
the plasma out of the disk via magnetic reconnections. 
The dusty torus and BLR in a RQQ are not shown. 

• m ~ 0.03: the thin disk is made of an inner and outer 
part; the inner disk extends to around ISCO. The two- 
phase ADAF is clumpy, so the jet produced via a mixture 
of the BP and BZ mechanisms is discontinuous. The dusty 
torus and BLR in a BL Lac are not shown. 

• m < 0.01: "Q" refers to the very low luminosity quiescent 
state that normally displays a hard power-law spectrum. 
The large scale magnetic fields rotating with the truncated 
thin disk and thick ADAF channel the continuous plasma 
to form coUimated and continuous jets via the BP mech- 
anism. There is neither dusty torus nor BLR in LLAGNs 
fT79l. 
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FIG. 40. Unification scheme of accretion-outfiow connections of accreting BHs for different ranges of accretion rates, corre- 
sponding to different spectral states or different types of astrophysical systems. From very low to very high accretion rates 
in units of the Eddington rate rii, these states are: hard/quiescent state, hard/low state, thermal-dominated soft state, steep 
power-law state, slim disk, and NDAF disk; different types of AGNs with similar accretion-outflow structures are also labeled 
for comparison. The ubiquitous dusty tori and BLRs are absent in all BHXBs and microquasars. The lack of dusty tori can 
be easily understood since there is no dust supply in them. It is then plausible that the lack of BLRs is the consequence of no 
dusty tori, suggesting that a BLR may be formed out of the evaporated inner dusty torus by the anisotropic radiation of the 
accretion disk in an AGN [T7^ . 



revise it inevitably. Of particular importance is the re- 
liable determinations of the mass and spin of accreting 
BHs. The BH mass allows us to determine the accretion 
rate in units of Eddington rate, i.e., m, which is a key 
parameter to the unification scheme of accretion-outflow 
connections. The BH spin is of course another key pa- 
rameter, because it can determine the radiative efhciency 
of the disk and jet power in the BZ mechanism. 



Once we have well determined BH parameters and a 
consistent and predictable theory describing the observed 
accretion-outflow connections, we are then ready to study 
the fundamental physics of BHs with BHXBs and micro- 
quasars, such as the properties of event horizon, space- 
time around Schwarzschild and Kerr BHs, BH spin en- 
ergy (and mass) extraction. For example, both the broad 
iron line and CF fitting methods can be used to measure 
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a BH's spin. However the two methods are equivalent 
only if the metric is accurately described with Kerr met- 
ric. Different spin measurements for a BH with the two 
different methods would invalidate Kerr metric, provided 
that we understand the accretion disk physics thoroughly. 
On the other hand, a correct metric is of course required 
to describe accurately the accretion flow around BHs. 
This is one example of the entanglement or interplay be- 
tween astrophysics and fundamental physics of accreting 
BHs. The stake is high, but the job is difficult. 



[1] R. Fender and T. Belloni. Stellar-Mass Black Holes and 
Ultraluminous X-ray Sources. Science, 337(6094):540- 
544, August 2012. 

[2] Ronald a. Romillard and Jeffrey E. McClintock. X-Ray 
Properties of Black-Hole Binaries. Annual Review of 
Astronomy and Astrophysics, 44(l):49-92, September 
2006. 

[3] Chris Done. Observational characteristics of accretion 

onto black holes. To appear in XXI Canary Islands 
Winter School of Astrophysics. Ed T Shahbaz, CUP, 
page 56, August 2010. 
[4] Chris Done, Marek Gierliski, and Aya Kubota. Mod- 
elling the behaviour of accretion flows in X-ray binaries. 
The Astronomy and Astrophysics Review, 15(l):l-66, 
August 2007. 

[5] Shuang-nan Zhang. Astrophysical Black Holes in the 
Physical Universe. In Zhang S. N. York, D. G., Gin- 
gerich, O., editor. Astronomy Revolution: 400 Years 
of Exploring the Cosmos, chapter Astrophysi, page 23. 
Taylor & Francis Group LLC/CRC Press, March 2011. 

[6] Ya B Zcl'dovich and Igor D Novikov. RELATIVISTIC 
ASTROPHYSICS. II. Soviet Physics Uspekhi, 8(4):522- 
577, April 1966. 

[7] Ya B Zel'dovich and Igor D Novikov. RELATIVISTIC 
ASTROPHYSICS. I. Soviet Physics Uspekhi, 7(6):763- 
788, June 1965. 

[8] B. LOUISE WEBSTER and PAUL MURDIN. Cygnus 
X-la Spectroscopic Binary with a Heavy Companion ? 
Nature, 235(5332) :37-38, January 1972. 

[9] I. F. Mirabel, L. F. Rodriguez, B. Cordicr, J. Paul, 
and F. Lebrun. A double-sided radio jet from the com- 
pact Galactic Centre annihilator 1E1740. 72942. Nature, 
358(6383) :215-217, July 1992. 
[10] I. F. Mirabel and L. F. Rodriguez. A superluminal 
source in the Galaxy. Nature, 371 (6492) :46-48, Septem- 
ber 1994. 

[11] S N Zhang, C A Wilson, B A Harmon, G J Fishman, 
W S Paciesas, and M Scott. lAUC 6046 : X-RAY N IN 
Sco ; 1994W. lAU Circ, 6046(1), 1994. 

[12] S. J. Tingay, D. L. Jauneey, R. A. Preston, J. E. 
Reynolds, D. L. Meier, D. W. Murphy, A. K. Tzioumis, 
D. J. McKay, M. J. Kesteven, J. E. J. LovcU, 
D. Campbell-Wilson, S. P. ElUngsen, R. Gough, R. W. 
Hunstead, D. L. Jonos, P. M. McCuIIoch, V. Migenes, 
J. Quick, M. W. Sinclair, and D. Smits. Relativis- 
tic motion in a nearby bright X-ray source. Nature, 
374(6518):141-143, March 1995. 



Acknowledgment. I appreciate inputs from Profs. Li- 
jun Gou, Wcimin Gu, Lixin Li, Bifang Liu, Dingxiong 
Wang, Wcnfei Yu, Feng Yuan, and Shu Zhang. The ed- 
itors of this book are thanked for inviting me to write 
this article, as well as their patience, persistency, proof 
reading it, and offering comments and suggestions to im- 
prove it in the end. This work is partially supported 
with funding by the 973 Program of China under grant 
2009CB824800, by the National Natural Science Founda- 
tion of China under grant Nos. 11133002 and 10725313, 
and by the Qianren start-up grant 292012312D1117210. 



[13] B. A. Harmon, C. A. Wilson, S. N. Zhang, W. S. Pa- 
ciesas, G. J. Fishman, R. M. Hjellming, M. P. Rupen, 
D. M. Scott, M. S. Briggs, andB. C. Rubin. Correlations 

between X-ray outbursts and relativistic ejections in the 
X-ray transient GRO J1655 40. Nature, 374(6524) :703- 
706, April 1995. 
[14] I F Mirabel and L. F. Rodriguez. SOURCES OF R EL- 
ATIVISTIC J ETS IN THE G ALAXY. Annual Review 
of Astronomy and Astrophysics, 37(l):409-443, Septem- 
ber 1999. 

[15] Rob Fender, Kinwali Wu, Helen Johnston, Tasso 
Tzioumis, Peter Jonker, Ralph Spencer, and Michiel 
Van Der Klis. An ultra-relativistic outflow from a neu- 
tron star accreting gas from a companion. Nature, 
427(6971) :222-4, January 2004. 

[16] J. F. Hao and S. N. Zhang. Large-Scale Cavities 
Surrounding Microquasars Inferred From Evolution of 
Their Relativistic Jets. The Astrophysical Journal, 
702(2):1648-1661, September 2009. 

[17] D. A. Kniffen, W. C. K. Alberts, D. L. Bertsch, B. L. 
Dingus, J. A. Esposito, C. E. Fichtel, R. S. Foster, 
R. C. Hartmau, S. D. Hunter, G. Kaiibacli, Y. C. Lin, 
J. R. Mattox, H. A. MayerHasscl wander, P. F. Michel- 
son, C. von Montigny, R. Mukherjee, P. L. Nolan, J. M. 
Paredes, P. S. Ray, E. J. Schneid, P. Sreekumar, M. Ta- 
vani, and D. J. Thompson. EGRET Observations of 
the GammaRay Source 2CG 135-fOl. The Astrophysi- 
cal Journal, 486(1):126-131, September 1997. 

[18] Josep M Paredes, Josep Mart, and Maria Massi. Discov- 
ery of a High-Energy Gamma-Ray Emitting Persistent 
Microquasar. Science, 288(June):2340-2342, 2000. 

[19] F Aharonian, A G Akhperjanian, K-M Aye, A R Bazer- 
Bachi, M Beilicke, W Benbow, D Berge, P Berghaus, 
K Bernlohr, C Boisson, O Bolz, V Borrel, I Braun, 
F Breitling, A M Brown, J Bussons Gordo, P M Chad- 
wick, L-M Chounet, R Cornils, L Costamantc, B De- 
grange, H J Dickinson, A Djannati-Atai', L O'c Drury, 
G Dubus, D Emmanoulopoulos, P Espigat, F Feinstein, 
P Fleury G Fontaine, Y Fuchs, S Funk, Y A Gallant, 
B Giebels, S Gillesscn, J F Glicenstein, P Goret, C Had- 
jichristidis, M Hauser, G Heinzelmann, G Henri, G Her- 
mann, J A Hinton, W Hofmann, M Holleran, D Horns, 
A Jacholkowska, O C de Jager, B Khelifl, Nu Komin, 
A Konopeiko, I J Latham, R Le Gallou, A Lemiere, 
M Lemoine-Goumard, N Leroy, T Lohse, A Marcowith, 
J-M Martin, O Martineau-Huynh, C Masterson, T J L 
McComb, M de Naurois, S J Nolan, A Noutsos, K J 



28 



Orford, J L Osborne, M Ouchrif, M Pantcr, G Pel- 
letier, S Pita, G Piihlhofer, M Punch, B C Rauben- 
heimer, M Raue, J Raux, S M Rayncr, A Rcimcr, 
O Reimer, J Ripken, L Rob, L RoUand, G Row- 
ell, V Sahakian, L Sauge, S Schlenker, R Schlickeiser, 
C Schuster, U Schwanke, M Siewert, H Sol, D Spangler, 
R Steenkamp, C Stegmann, J-P Tavernet, R Terrier, 
C G Thcoret, M Tluczykont, G Vasilciadis, C Venter, 
P Vincent, H J Volk, and S J Wagner. Discovery of very 
high energy gamma rays associated with an x-ray bi- 
nary. Science (New York, N.Y.), 309(5735):746-9, July 
2005. 

[20] F Aharonian. A New Population of Very High En- 
ergy Gamma-Ray Sources in the Milky Way. Science, 
307(Marcli): 1938- 1942, 2005. 

[21] 1. F. Mirabel. Gamma-ray binaries. Astrophysics and 
Space Science, 309(l-4):267-270, May 2007. 

[22] Wei Cui. TeV gamma-ray astronomy. Research in As- 
tronomy and Astrophysics, 9(8):841-860, August 2009. 

[23] Gustavo E. Romero. Gamma-ray emission from 
pulsar/massive-star binaries. Proceedings of the Inter- 
national Astronomical Union, 5(H15):126-130, October 
2010. 

[24] Jian Li, Diego F. Torres, Shu Zhang, Daniela Hadasch, 
Nanda Rea, G. Andrea Caliandro, Yupeng Chen, and 
Jianmin Wang. Unveiling the Super-Orbital Modulation 

of Ls 1 -1-61303 in X-Rays. The Astrophysical Journal, 
744(1):L13, .Jarmary 2012. 

[25] K. S. Thornc. Classical Black Holes: The Nonlinear Dy- 
namics of Curved Spacetime. Science, 337(6094):536- 
538, August 2012. 

[26] E. Wittcn. Quantum Mechanics of Black Holes. Science, 
337(6094) :538-540, August 2012. 

[27] M. Volontcri. The Formation and Evolution of Massive 
Black Holes. Science, 337(6094):544-547, August 2012. 

[28] Tomaso M Belloni, Sara E Motta, and Teodoro Munoz 
Darias. Black hole transients. Bulletin of the Astronmi- 
cal Society of India, 39(3):409-428, September 2011. 

[29] J. N. Zhou, Q. Z. Liu, Y. P. Chen, J. Li, J. L. Qu, 
S. Zhang, H. Q. Gao, and Z. Zhang. The last three 
outbursts of H1743-322 observed by RXTE in its latest 
service phase. Monthly Notice of the Royal Astronomical 
Society, 2013. 

[30] Andrew G. CantrcU, Charles D. Bailyn, Jerome a. 
Orosz, Jeffrey E. McClintock, Ronald a. Rcmillard, 
Cynthia S. Froning, Joseph Neilsen, Dawn M. Gelino, 
and Lijun Gou. the Inclination of the Soft X-Ray Tran- 
sient a062000 and the Mass of Its Black Hole. The As- 
trophysical Journal, 710(2):1127-1141, February 2010. 

[31] Lijun Gou, Jeffrey E. McClintock, James F. Steiner, 
Ramcsh Narayan, Andrew G. Cantroll, Charles D. Bai- 
lyn, and Jerome a. Orosz. the Spin of the Black Hole in 
the Soft X-Ray Transient a0620-00. The Astrophysical 
Journal, 718(2) :L122-L126, August 2010. 

[32] Rebecca Shafee, Jeffrey E. McClintock, Ramesh 
Narayan, Shane W Davis, Li-xin Li, and Ronald A 
Remillard. Estimating the Spin of Stellar-Mass Black 
Holes by Spectral Fitting of the X-Ray Continuum. 
The Astrophysical Journal, 636(2):L113-L116, January 
2006. 

[33] James F. Steiner, Rubens C. Reis, Jeffrey E. McClin- 
tock, Ramesh Narayan, Ronald a. Remillard, Jerome a. 
Orosz, Lijun Gou, Andrew C. Fabian, and Manuel a. P. 
Torres. The spin of the black hole microquasar XTE 



J1550564 via the continuum-fitting and Fo-line meth- 
ods. Monthly Notices of the Royal Astronomical Society, 
416(2):941 958, September 2011. 

[34] James F. Steiner, Jeffrey E. McChntock, and Mark J. 
Reid. THE DISTANCE, INCLINATION, AND SPIN 
OF THE BLACK HOLE MICROQUASAR H1743-322. 
The Astrophysical Journal, 745(1):L7, January 2012. 

[35] Jeffrey E. McClintock, Rebecca Shafee, Ramesh 
Narayan, Ronald A. Remillard, Shane W. Davis, and 
LiXin Li. The Spin of the NearExtrcme Kerr Black Hole 
GRS 1915-M05. The Astrophysical Journal, 652(1):518- 
539, November 2006. 

[36] Jerome a. Orosz, Jeffrey E. McClintock, Jason P. Auf- 
denberg, Ronald a. Remillard, Mark J. Reid, Ramesh 
Narayan, and Lijun Gou. the Mass of the Black Hole 
in Cygrms X-1. The Astrophysical Journal, 742(2) :84, 
December 2011. 

[37] Lijun Gou, Jeffrey E. McClintock, Mark J. Reid, 
Jerome a. Orosz, James F. Steiner, Ramesh Narayan, 
Jingen Xiang, Ronald a. Remillard, Keith a. Ar- 
naud, and Shane W. Davis, the Extreme Spin of the 
Black Hole in Cygnus X-1. The Astrophysical Journal, 
742(2):85, December 2011. 

[38] Shane W. Davis, Chris Done, and Omer M. Blaes. Test- 
ing Accretion Disk Theory in Black Hole XRay Bina- 
ries. The Astrophysical Journal, 647(1) :525-538, August 
2006. 

[39] Jerome a. Orosz, Danny Steeghs, Jeffrey E. McClin- 
tock, Manuel a. P. Torres, Ivan Bochkov, Lijun Gou, 
Ramesh Narayan, Michael Blaschak, Alan M. Levine, 
Ronald a. Remillard, Charles D. Bailyn, Morgan M. 
Dwyer, and Michelle Buxton, a New Dynamical Model 
for the Black Hole Binary Lmc X-1. The Astrophysical 
Journal, 697(1):573 591, May 2009. 

[40] Lijun Gou, Jeffrey E. McClintock, Jifeng Liu, Ramesh 
Narayan, James F. Steiner, Ronald a. Remillard, 
Jerome a. Orosz, Shane W. Davis, Ken Ebisawa, and 
Eric M. Schlegel. a Determination of the Spin of the 
Black Hole Primary in Lmc X-1. The Astrophysical 
Journal, 701(2):1076-1090, August 2009. 

[41] James F Steiner, Rubens C Reis, Andrew C Fabian, 
Ronald A Remillard, Jeffrey E. McClintock, Lijun Gou, 
Ryan Cooke, Laura W Brenneman, and Jeremy S 
Sanders. A broad iron line in LMC X-1. Monthly No- 
tices of the Royal Astronomical Society, 427(3):2552- 
2561, December 2012. 

[42] A H Prestwich, R Kilgard, P A Crowther, S Carpano, 
A M T Pollock, A Zezas, S H Saar, T P Roberts, and 
M J Ward. The Orbital Period of the Wolf-Rayet Binary 
IC 10 X-1: Dynamic Evidence that the Compact Object 
Is a Black Hole. The Astrophysical Journal, 669(1) :L21- 
L24, November 2007. 

[43] Jeffrey M Silverman and Alexei V Filippenko. On IC 10 
X-1, the Most Massive Known Stellar-Mass Black Hole. 
The Astrophysical Journal, 678(1):L17-L20, May 2008. 

[44] P. a. Crowther, R. Barnard, S. Carpano, J. S. Clark, 
V. S. Dhillon, and a. M. T. Pollock. NGC 300 X-1 is a 
Wolf-Rayet/black hole binary. Monthly Notices of the 
Royal Astronomical Society: Letters, 403(1):L41-L45, 
March 2010. 

[45] Jerome a Orosz, Jeffrey E McClintock, Ramesh 
Narayan, Charles D Bailyn, Joel D Hartman, Lucas 

Macri, Jiefeng Liu, Wolfgang Pictsch, Ronald a Remil- 
lard, Avi Shporer, and Tsevi Mazeh. A 15.65-solar-mass 



29 



black hole in an eclipsing binary in the nearby spiral 
galaxy M 33. Nature, 449(7164) :872-5, October 2007. 

[46] Jifeng Liu, Jeffrey E. McClintock, Ramesh Narayan, 
Shane W. Davis, and Jerome a. Orosz. ERRATUM: 
PRECISE MEASUREMENT OF THE SPIN PARAM- 
ETER OF THE STELLAR-MASS BLACK HOLE M33 
X-7 (2008, ApJ, 679, L37). The Astrophysical Journal, 
719(1):L109-L109, August 2010. 

[47] S. N. Zhang, Wei Cui, and Wan Chen. Black Hole Spin 
in X-Ray Binaries: Observational Consequences. The 
Astrophysical Journal, 482(2) :L155-L158, June 1997. 

[48] Jeffrey E McClintock, Ramesh Narayan, Shane W 
Davis, Lijun Gou, Akshay Kulkarni, Jerome A Orosz, 
Robert F Penna, Ronald A Remillard, and James F 
Steiner. Measuring the spins of accreting black holes. 
Classical and Quantum Gravity, 28(11):114009, June 
2011. 

[49] Jean-pierre Lasota. The disc instability model of dwarf 
novae and low-mass X-ray binary transients. New As- 
tronomy Reviews, 45:1993-1995, 2001. 

[50] Jean-Pierre Lasota. Disc outbursts in various types of 
binary systems. Memorie della Societa Astronomica 
Italiana, 83:469, November 2012. 

[51] N. I. Shakura and R A Sunyaev. Black Holes in Bi- 
nary Systems. Observational Appearance. Astronomy 
and Astrophysics, 24:337-355, 1973. 

[52] Don N. Page and Kip S. Thorne. Disk- Accretion onto a 
Black Hole. Time- Averaged Structure of Accretion Disk. 
The Astrophysical Journal, 191:499, July 1974. 

[53] Joseph Neilsen and Julia C Lee. Accretion disk winds as 
the jet suppression mechanism in the microquasar GRS 
1915-M05. Nature, 458(7237) :481-4, March 2009. 

[54] G Ponti, R P Fender, M C Begelman, R J H Dunn, 
J Neilsen, and M Coriat. Ubiquitous equatorial accre- 
tion disc winds in black hole soft states. Monthly Notices 
of the Royal Astronomical Society: Letters, 422(1):L11- 
L15, May 2012. 

[55] M. A. Abramowicz and W. Kluzniak. A precise determi- 
nation of black hole spin in GRO J1655-40. Astronomy 
and Astrophysics, 374(3) :L19-L20, August 2001. 

[56] G. Torok, M. A. Abramowicz, W. Kluzniak, and 
Z. Stuchlfk. The orbital resonance model for twin peak 
kHz quasi periodic oscillations in microquasars. Astron- 
omy and Astrophysics, 436(1):1 8. June 2005. 

[57] Chang- Yin Huang, Zhao-Ming Gan, Jiu-Zhou Wang, 
and Ding-Xiong Wang. A resonance model with mag- 
netic connection for 3:2 HFQPO pairs in black hole bina- 
ries. Monthly Notices of the Royal Astronomical Society, 
403(4): 1978-1982, April 2010. 

[58] Mari Kolehmainen and Chris Done. Limits on spin 
determination from disc spectral fitting in GX 3394. 
Monthly Notices of the Royal Astronomical Society, 
406(4):2206-2212, August 2010. 

[59] Ramesh Narayan and Insu Yi. Advection-dominated 
Accretion: Underfed Black Holes and Neutron Stars. 
The Astrophysical Journal, 452:710, October 1995. 

[60] B. F. Liu, W. Yuan, F. Meyer, E. Meyer-Hofmeister, 
and G. Z. Xie. Evaporation of Accretion Disks around 
Black Holes: The Disk-Corona Transition and the Con- 
nection to the Advection-dominated Accretion Flow. 
The Astrophysical Journal, 527(1) :L17-L20, December 
1999. 

[61] J.M. Miller. Relativistic X-Ray Lines from the Inner 
Accretion Disks Around Black Holes. Annual Review of 



Astronomy and Astrophysics, 45(l):441-479, September 
2007. 

[62] A R King, U Kolb, and E Szuszkiewicz. Why LowMass 
Black Hole Binaries Are Transient. The Astrophysical 
Journal, 488(l):89-93, October 1997. 

[63] D. Altamirano, T. Belloni, M. Linares, M. van der Klis, 
R. Wijnands, P. a. Curran, M. Kalamkar, H. Stiele, 
S. Motta, T. Munoz Darias, P. Casella, and H. Krimm. 
the Faint Heartbeats of Igr J170913624: an Excep- 
tional Black Hole Candidate. The Astrophysical Jour- 
nal, 742(2):L17, December 2011. 

[64] Shan-shan Weng and Shuang-nan Zhang. COM- 
PARING THE ACCRETION DISK EVOLUTION OF 
BLACK HOLE AND NEUTRON STAR X-RAY BL 
NARIES FROM LOW TO SUPER-EDDINGTON LU- 
MINOSITY. The Astrophysical Journal, 739(1):42, 
September 2011. 

[65] Yangsen Yao, S. Nan Zhang, Xiaoling Zhang, Yuxin 
Feng, and Craig R. Robinson. Studying the Properties 
of Accretion Disks and Coronae in Black Hole XRay Bi- 
naries with Monte Carlo Simulation. The Astrophysical 
Journal, 619(l):446-454, January 2005. 

[66] Shuang Nan Zhang. Similar phenomena at different 
scales: black holes, the Sun, 7-ray bursts, supernovae, 
galaxies and galaxy clusters. Highlights of Astronomy, 
14:14-25, December 2007. 

[67] Lin Wang, Cheng Fang, and Yu Ying Liu. Flux-rms 
relation in solar radio bursts. Astrophysics and Space 
Science, 318(l-2):79-86, September 2008. 

[68] J Biteau and B Giebels. The minijets-in-a-jet statisti- 
cal model and the rms-flux correlation. Astronomy & 
Astrophysics, 548:A123, December 2012. 

[69] J M Bardeen. Kerr metric black holes. Nature, 
226(5240):64 5, April 1970. 

[70] Kip S. Thorne. Disk-Accretion onto a Black Hole. 
II. Evolution of the Hole. The Astrophysical Journal, 
191:507, July 1974. 

[71] a. R. King and U. Kolb. The evolution of black hole 
mass and angular momentum. Monthly Notices of the 
Royal Astronomical Society, 305(3) :654 660, May 1999. 

[72] Ramesh Narayan and Jeffrey E. McClintock. Observa- 
tional evidence for a correlation between jet power and 
black hole spin. Monthly Notices of the Royal Astronom- 
ical Society: Letters, 419(1):L69 L73, January 2012. 

[73] James F. Steiner, Jeffrey E. McClintock, and Ramesh 
Narayan. Jet Power and Black Hole Spin: Testing an 
Empirical Relationship and Using It To Predict the 
Spins of Six Black Holes. The Astrophysical Journal, 
762(2):104, January 2013. 

[74] R. P. Fender, E. Gallo, and D. Russell. No evidence 
for black hole spin powering of jets in X-ray bina- 
ries. Monthly Notices of the Royal Astronomical Society, 
1434:no-no, July 2010. 

[75] R. D. Blandford and D. G. Payne. Hydromagnetic 
flows from accretion discs and the production of radio 
jets. Monthly Notices of the Royal Astronomical Society, 
199:883-903, 1982. 

[76] R. D. Blandford and R. L. Znajek. Electromagnetic 
extraction of energy from Kerr black holes. Monthly 
Notices of the Royal Astronomical Society, 179:433-456, 
1977. 

[77] Li Xue, Aleksander Sdowski, Marek a. Abramowicz, and 
Ju-Fu Lu. Studies of Thermally Unstable Accretion 
Disks Around Black Holes With Adaptive Pseudospec- 



30 



tral Domain Decomposition Method. li. Limit-Cycle Be- 
havior in Accretion Disks Around Kerr Black Holes. The 
Astrophysical Journal Supplement Series, 195(1):7, July 
2011. 

[78] James M. Bardeen, William H. Press, and Saul A. 
Teukolsky. Rotating Black Holes: Locally Nonrotating 
Frames, Energy Extraction, and Scalar Synchrotron Ra- 
diation. The Astrophysical Journal, 178:347, December 
1972. 

[79] S N Zhang, K Ebisawa, R Sunyaev, Y Ueda, B A Har- 
mon, S Sazonov, G J Fishman, H Inoue, W S Paciesas, 
and T Takahash. Broadband HighEnergy Observations 

of the Suporluminal Jet Source GRO J 165540 during an 
Outburst. The Astrophysical Journal, 479(1):381 387, 
April 1997. 

[80] Jerome A Orosz and Charles D Bailyn. Optical Ob- 
servations of GRO J165540 in Quiescence. I. A Precise 
Mass for the Black Hole Primary. The Astrophysical 
Journal, 477(2) :876-896, March 1997. 

[81] LiXin Li, Erik R. Zimmerman, Ramcsh Narayan, and 
Jeffrey E. McClintock. Multitemperature Blackbody 
Spectrum of a Thin Accretion Disk around a Kerr Black 
Hole: Model Computations and Comparison with Ob- 
servations. The Astrophysical Journal Supplement Se- 
mes, 157(2):335-370, April 2005. 

[82] Shane W Davis and Ivan Hubeny. A Grid of Relativis- 
tic, NonLTE Accretion Disk Models for Spectral Fitting 
of Black Hole Binaries. The Astrophysical Journal Sup- 
plement Series, 164(2):530 535, June 2006. 

[83] Jeffrey E. McChntock and Ronald A Remillard. Mea- 
suring the Spins of Stellar-Mass Black Holes. Astro2010 
Science White Paper, February 2009. 

[84] Robert F. Pcnna, Jonathan C. McKinncy, Ramesh 
Narayan, Alexander Tctickhovskoy, Rebecca Shafee, 
and Jeffrey E. McClintock. Simulations of magnetized 
discs around black holes: effects of black hole spin, disc 
thickness and magnetic field geometry. Monthly Notices 
of the Royal Astronomical Society, 408(2):752-782, Oc- 
tober 2010. 

[85] Akshay K. Kulkarni, Robert F. Penna, Roman V. 
Shcherbakov, James F. Steiner, Ramesh Narayan, Alek- 
sander Sdowski, Yucong Zhu, Jeffrey E. McClintock, 
Shane W. Davis, and Jonathan C. McKinney. Measur- 
ing black hole spin by the continuum-fitting method: ef- 
fect of deviations from the Novikov-Thorne disc model. 
Monthly Notices of the Royal Astronomical Society, 
414(2):1183-1194, June 2011. 

[86] Yucong Zhu, Shane W. Davis, Ramesh Narayan, Ak- 
shay K. Kulkarni, Robert F. Penna, and Jeffrey E. Mc- 
Clintock. The eye of the storm: light from the inner 
plunging region of black hole accretion discs. Monthly 
Notices of the Royal Astronomical Society, 424(4):2504- 
2521, August 2012. 

[87] James F. Steiner, Jeffrey E. McChntock, Ronald a. 
Remillard, Ramesh Narayan, and Lijun Gou. MEA- 
SURING BLACK HOLE SPIN VIA THE X-RAY 
CONTINUUM-FITTING METHOD: BEYOND THE 
THERMAL DOMINANT STATE. The Astrophysical 
Journal, 701(2):L83-L86, August 2009. 

[88] James F. Steiner, Jeffrey E. McClintock, Ronald a. 
Remillard, Lijun Gou, Shin'ya Yamada, and Ramesh 
Narayan. the Constant Inner-Disk Radius of Lmc X-3: 
a Basis for Measuring Black Hole Spin. The Astrophys- 
ical Journal, 718(2):L117-L121, August 2010. 



[89] Dacheng Lin, Ronald a. Remillard, and Jeroen Homan. 
Spectral States of Xte J1701 462: Link Between Z and 
Atoll Sources. The Astrophysical Journal, 696(2):1257- 
1277, May 2009. 

[90] Shan-shan Weng, Shuang-nan Zhang, and Hai-hui Zhao. 
Super-Eddington Accretion in the Ultraluminous X-ray 
Source NGC1313 X-2: An Ephemeral Feast. The As- 
trophysical Journal, submitted, 2013. 

[91] Wei-Min Gu. Radiation Pressure-Supported Accretion 
Disks: Vertical Structure, Energy Advcction, and Con- 
vective Stability. The Astrophysical Journal, 753(2): 118, 
July 2012. 

[92] Xie Chen, Shuang Nan Zhang, and Guo Qiang Ding. 
PropcUerdriven Spectral State Transition in the Low- 
Mass XRay Binary 4U 160852. The Astrophysical Jour- 
nal, 650(l):299-302, October 2006. 

[93] Feng Yuan, Wei Cui, and Ramesh Narayan. An Accre- 
tionjet Model for Black Hole Binaries: Interpreting the 
Spectral and Timing Features of XTE J1118-I-480. The 
Astrophysical Journal, 620(2) :905-914, February 2005. 

[94] S. N. Zhang. The 1996 soft state transition of cygnus 
x-1. ApJ, 477:L95-98, 1997. 

[95] Qingjuan Yu and Scott Tremaine. Observational con- 
straints on growth of massive black holes. Monthly No- 
tices of the Royal Astronomical Society, 335(4):965-976, 
October 2002. 

[96] Tsing-Wai Wong, Franccsca Valsecchi, Tassos Fragos, 
and Vassiliki Kalogcra. UNDERSTANDING COM- 
PACT OBJECT FORMATION AND NATAL KICKS. 
HI. THE CASE OF CYGNUS X-1. The Astrophysical 
Journal, 747(2):111, March 2012. 
[97] James F Steiner, Ramesh Narayan, Jeffrey E. McClin- 
tock, and Ken Ebisawa. A Simple Comptonization 
Model. Publications of the Astronomical Society of the 
Pacific, 121(885):1279-1290, November 2009. 
[98] Andrew G Cantrell, Charles D Bailyn, Jeffrey E. Mc- 
Clintock, and Jerome A Orosz. Optical State Changes 
in the X-Ray-quiescent Black Hole A0620-00. The As- 
trophysical Journal, 673(2):L159-L162, February 2008. 
[99] Laura Kreidberg, Charles D. Bailyn, Will M. Farr, and 
Vicky Kalogera. Mass Measurements of Black Holes in 
X-Ray Transients: Is There a Mass Gap? The Astro- 
physical Journal, 757(1) :36, September 2012. 

[100] F Van Dcr Hooft, F Alberts, and J Van Paradijs. AS- 
TRONOMY AND The quiescence optical light curve 
of Nova Scorpii 1994. Astronomy and Astrophysics, 
550:538-550, 1998. 

[101] S N Philhps, T Shahbaz, and Ph Podsiadlowski. The 
outburst radial velocity curve of X-ray Nova Scorpii 
1994 (= GRO J1655-40): a reduced mass for the black 
hole? Monthly Notices of the Royal Astronomical Soci- 
ety, 304(4) :839 844, April 1999. 

[102] T Shahbaz, F Van Der Hooft, J Casares, P A Charles, 
and J Van Paradijs. The mass of X-ray Nova Scorpii 
1994 (=GRO J1655-40). Monthly Notices of the Royal 
Astronomical Society, 306(1) :89-94, June 1999. 

[103] Jenny Greene, Charles D Bailyn, and Jerome A Orosz. 
Optical and Infrared Photometry of the Microquasar 
GRO J165540 in Quiescence. The Astrophysical Jour- 
nal, 554(2):1290-1297, June 2001. 

[104] Martin E. Beer and Philipp Podsiadlowski. The qui- 
escent light curve and the evolutionary state of GRO 
J1655-40. Monthly Notices of the Royal Astronomical 
Society, 331(2):351-360, March 2002. 



31 



[105] T. Shahbaz. Determining the spectroscopic mass ratio 
in interacting binaries: application to X-Ray Nova Sco 
1994. Monthly Notices of the Royal Astronomical Soci- 
ety, 339(4): 1031-1040, March 2003. 

[106] James M. Bardeen and Jacobus A. Petterson. The 
Lense-Thirring Effect and Accretion Disks around Kerr 
Black Holes. The Astrophysical Journal, 195:L65, Jan- 
uary 1975. 

[107] R. M. Hjellming and M. P. Rupon. Episodic ejection of 
relativistic jets by the X-ray transient GRO J1655 - 40. 
Nature, 375(6531) :464-468, June 1995. 

[108] C Foellmi, E Depagne, T H Dall, and I F Mirabel. On 
the distance of GRO J1655-40. Astronomy and Astro- 
physics, 457(l):249-255, October 2006. 

[109] Ccdric Foellmi. What is the closest black hole to the 
Sun? New Astronomy, 14(8):674-691, November 2009. 

[110] D. Steeghs and J. Casares. The Mass Donor of Scorpius 
XI Revealed. The Astrophysical Journal, 568(1):273- 
278, March 2002. 

[Ill] Carole A. Haswell and Allen W. Shafter. A detection 
of orbital radial velocity variations of the primary com- 
ponent in the black hole binary A0620 - 00 (= V616 
Monocerotis). The Astrophysical Journal, 359:L47, Au- 
gust 1990. 

[112] Jerome A. Orosz, Charles D. Bailyn, Ronald A. Remil- 
lard, Jeffrey E. McClintock, and Craig B. Foltz. Qui- 
escent accretion disks in black hole X-ray novae. The 
Astrophysical Journal, 436:848, December 1994. 

[113] Roberto Soria, Dayal T. Wickramasinghe, Richard W. 
Hunstead, and Kinwah Wu. Measuring the Motion of 
the Black Hole in GRO J165540. The Astrophysical 
Journal, 495(2):L95-L98, March 1998. 

[114] Shuang-Nan Zhang, Jinyuan Liao, and Yangsen Yao. 
Measuring the black hole masses in accreting X-ray bi- 
naries by detecting the Doppler orbital motion of their 
accretion disc wind absorption lines. Monthly Notices 
of the Royal Astronomical Society, 421(4):3550-3556, 
April 2012. 

[115] Jon M Miller, John Raymond, Andy Fabian, Danny 
Steeghs, Jeroen Homan, Chris Reynolds, Michiel 
van der Klis, and Rudy Wijnands. The magnetic 
nature of disk accretion onto black holes. Nature, 
441 (7096) :953-5, June 2006. 

[116] P. A. CONNORS and R. F. STARK. Observable gravi- 
tational effects on polarised radiation coming from near 
a black hole. Nature, 269(5624):128-129, September 
1977. 

[117] P. A. Connors, R. F. Stark, and T. Piran. Polarization 
features of X-ray radiation emitted near black holes. 
The Astrophysical Journal, 235:224, January 1980. 

[118] Li-Xin Li, Ramosh Narayan, and Jeffrey E. McClintock. 
Inferring the Inclination of a Black Hole Accretion Disk 
From Observations of Its Polarized Continuum Radia- 
tion. The Astrophysical Journal, 691(l):847-865, Jan- 
uary 2009. 

[119] Jeremy D. Schnittman and Julian H. Krolik. X-Ray 
Polarization From Accreting Black Holes: the Thermal 
State. The Astrophysical Journal, 701(2):1175-1187, 
August 2009. 

[120] Subrahmanyan Chandrasekhar. Radiative transfer. 
Dover, New York, 1960. 

[121] J M Miller, C S Reynolds, A C Fabian, E M Cack- 
ett, G Miniutti, J Raymond, D Steeghs, R Reis, and 
J Homan. Initial Measurements of Black Hole Spin in 



GX 339-4 from Suzaku Spectroscopy. The Astrophysical 
Journal, 679(2):L113-L116, June 2008. 
[122] J. M. Miller, C. S. Reynolds, a. C. Fabian, G. Miniutti, 
and L. C. Gallo. Stellar-Mass Black Hole Spin Con- 
straints From Disk Reflection and Continuum Model- 
ing. The Astrophysical Journal, 697(1):900-912, May 
2009. 

[123] R. C. Reis, a. C. Fabian, R. R. Ross, and J. M. Miller. 
Determining the spin of two stellar-mass black holes 
from disc reflection signatures. Monthly Notices of 
the Royal Astronomical Society, 395(3): 1257-1264, May 
2009. 

[124] Zhixing Ling, Shuang Nan Zhang, Jingen Xiang, and 
Shichao Tang, a Study of the X-Ray Dust Scattering 
Halo of Cyg X-1 With a Cross-Correlation Method. The 
Astrophysical Journal, 690(l):224-230, January 2009. 

[125] Zhixing Ling, Shuang Nan Zhang, and Shichao Tang. 
DETERMINING THE DISTANCE OF Cyg X-3 WITH 
ITS X-RAY DUST SCATTERING HALO. The Astro- 
physical Journal, 695(2):1111-1120, April 2009. 

[126] Mark J. Reid, Jeffrey E. McClintock, Ramesh Narayan, 
Lijun Gou, Ronald a. Remillard, and Jerome a. Orosz. 
the Trigonometric Parallax of Cygnus X-1. The Astro- 
physical Journal, 742(2):83, December 2011. 

[127] Bozena Czerny, Krzysztof Hryniewicz, M. Nikoajuk, 
and A. Sadowski. Constraints on the black hole spin 
in the quasar SDSS J094533. 99+100950.1. Monthly No- 
tices of the Royal Astronomical Society, pages no-no, 
June 2011. 

[128] Makoto Kishimoto, Robert Antonucci, Omer Blaes, 
Andy Lawrence, Catherine Boisson, Marcus Albrecht, 
and Christian Leipski. The characteristic blue spectra of 
accretion disks in quasars as uncovered in the infrared. 
Nature, 454(7203):492-4, July 2008. 

[129] Renyu Hu and Shuang-nan Zhang. Optical polariza- 
tion of quasars and the Balmer edge feature revealed by 
ultraviolet and polarized visible to near-infrared emis- 
sions. Monthly Notices of the Royal Astronomical Soci- 
ety, 426(4) :2847 2858, November 2012. 

[130] Hui Zhang, Feng Yuan, and Sylvain Chaty. Modeling 
the Hard States of Three Black Hole Candidates. The 
Astrophysical Journal, 717(2) :929-936, July 2010. 

[131] Feng Yuan, Zhaolong Yu, and Luis C. Ho. Revisiting 
the Fundamental Plane of Black Hole Activity At Ex- 
tremely Low Luminosities. The Astrophysical Journal, 
703(1):1034-1043, September 2009. 

[132] Feng Yuan and Wei Cui. RadioXRay Correlation and 
the Quiescent State of Black Hole Sources. The Astro- 
physical Journal, 629(1):408-413, August 2005. 

[133] Feng Yuan. Luminous Hot Accretion Flows: Thermal 
Equilibrium Curve and Thermal Stability. The Astro- 
physical Journal, 594(2) :L99-L102, September 2003. 

[134] Feng Yuan, Andrzej A Zdziarski, Yongquan Xue, and 
Xue-bing Wu. ModeUng the Hard States of XTE 
J1550564 during Its 2000 Outburst. The Astrophysical 
Journal, 659(1) :541-548, April 2007. 

[135] F Meyer and B F Liu. Evaporation: The change from 
accretion via a thin disk to a coronal flow. Astronomy 
and Astrophysics, 361:175-188, 2000. 

[136] F Meyer and B F Liu. Black hole X-ray binaries : a new 
view on soft-hard spectral transitions ASTRONOMY 
AND. Astronomy and Astrophysics, 354:L67-L70, 2000. 

[137] B. F. Liu, S. Minoshige, F. Meyer, E. MeyerHofmeis- 
ter, and T. Kawaguchi. TwoTemperature Coronal 



32 



Flow above a Thin Disk. The Astrophysical Journal, 
575(1):117-126, August 2002. 

[138] B F Liu, Ronald E Taam, E. McyorHofmcistcr, and 
F. Meyer. The Existence of Inner Cool Disks in the 
Low/Hard State of Accreting Black Holes. The Astro- 
physical Journal, 671(l):695-705, December 2007. 

[139] B. F. Liu, C. Done, and Ronald E. Taam. the Effect 
of Coronal Radiation on a Residual Inner Disk in the 
Low/Hard Spectral State of Black Hole X-Ray Binary 
Systems. The Astrophysical Journal, 726(1):10, January 
2011. 

[140] B F Liu, F Meyer, and E. Meyer-Hofmeister. An inner 

disk below the ADAF: the intermediate spectral state 
of black hole accretion. Astronomy and Astrophystcs, 
454(1):L9-L12, July 2006. 

[141] F Meyer, B F Liu, and E. Meyer-Hofmeister. Re- 
condensation from an ADAF into an inner disk: the 
intermediate state of black hole accretion? Astronomy 
and Astrophysics, 463(l):l-9, February 2007. 

[142] Ronald E Taam, B. F. Liu, F Meyer, and E. Meyer- 
Hofmeister. On the Properties of Inner Cool Disks in 
the Hard State of Black Hole XRay Transient Systems. 
The Astrophysical Journal, 688(1) :527-536, November 
2008. 

[143] Erlin Qiao and Bi Fang Liu. Dependence of Spectral 
State Transition and Disk Truncation on Viscosity Pa- 
rameter . Publications of the Astronomical Society of 
Japan, 61(2):403-410, 2009. 

[144] Erlin Qiao and B. F. Liu. the Emission From an Inner 
Disk and a Corona in the Low and Intermediate Spectral 
States of Black Hole X-Ray Binaries. The Astrophysical 
Journal, 744(2):145, January 2012. 

[145] B F Liu, F Meyer, E. Mcycr-Hofmcister, B F Liu, and 
F Meyer. Spectral state transitions in low-mass X-ray 
binaries the effect of hard and soft irradiation. Astron- 
omy and Astrophysics, 442(2):555-562, November 2005. 

[146] E. Meyer-Hofmeister, B F Liu, and F Meyer. Hysteresis 
in spectral state transitions ? a challenge for theoreti- 
cal modeling. Astronomy and Astrophysics, 432(1):181- 
187, March 2005. 

[147] E. Meyer-Hofmcister, B. F. Liu, and F. Meyer. The 
hard to soft spectral transition in LMXBs-affected by 
recondensation of gas into an inner disk. Astronomy 
and Astrophysics, 508(1):329 337, December 2009. 

[148] Jian-Min Wang, Cheng Cheng, and Yan-Rong Li. 
Clumpy Accretion Onto Black Holes. I. Clumpy- 
Advection-Dominated Accretion Flow Structure and 
Radiation. The Astrophysical Journal, 748(2):147, April 
2012. 

[149] S. N. Zhang, B. A. Harmon, W. S. Paciesas, G. J. Fish- 
man, J. E. Grindlay, D. Barret, M. Tavani, P. Kaaret, 
P. Bloser, E. Ford, and L. Titarchuk. Low state hard 
X-ray outburst from the X-ray burster 4U 1608-522 ob- 
served by BATSE/CGRO. Astronomy and Astrophysics 
Supplement, 120:279-282, 1996. 

[150] Yu-Peng Chen, Shu Zhang, Shuang-Nan Zhang, Jian 
Li, and Jian-Min Wang. Type-I X-Ray Bursts Reveal a 
Fast Co-Evolving Behavior of the Corona in an X-Ray 
Binary. The Astrophysical Journal, 752(2) :L34, June 
2012. 

[151] S. N. Zhang. Three-Layered Atmospheric Structure in 
Accretion Disks Around Stellar-Mass Black Holes. Sci- 
ence, 287(5456):1239-1241, February 2000. 

[152] Feng Yuan, Jun Lin, Kinwah Wu, and Luis C. Ho. 



A magnetohydrodynamical model for the formation of 
episodic jets. Monthly Notices of the Royal Astronomi- 
cal Society, 395(4):2183-2188, June 2009. 

[153] B Czerny, M. Niko?ajuk, A. R?a?ska, A.-M. Dumont, 
Z Loska, and P. T. ?ycki. Universal spectral shape of 
high accretion rate AGN. Astronomy and Astrophysics, 
412(2) :317-329, December 2003. 

[154] Y. L. Ai, W. Yuan, H. Y. Zhou, T. G. Wang, and 
S. H. Zhang. X-Ray Properties of Narrow-Line Seyfert 
1 Galaxies With Very Small Breadline Widths. The 
Astrophysical Journal, 727(1):31, January 2011. 

[155] Chris Done, S. W. Davis, C. Jin, O. Blaes, and M. Ward. 
Intrinsic disc emission and the soft X-ray excess in active 
galactic nuclei. Monthly Notices of the Royal Astronom- 
ical Society, 420(3): 1848-1860, March 2012. 

[156] Sergei Nayakshin and Fulvio Melia. Magnetic Flares and 
the Observed tau T approximately 1 in Seyfert Galaxies. 
The Astrophysical Journal, 490(1) :L13-L16, November 
1997. 

[157] Sergei Nayakshin and James B. Dove. Accretion Disks 
with Coronae in Cygnus XI: The Role of a Transition 
Layer. The Astrophysical Journal, 560(2):885-891, Oc- 
tober 2001. 

[158] Wenfei Yu, Michiel van der KUs, and Rob Fender. The 
Correlation between Hard X-Ray Peak Flux and Soft X- 
Ray Peak Flux in the Outburst Rise of Low-Mass X-Ray 
Binaries. The Astrophysical Journal, 611(2):L121-L124, 
August 2004. 

[159] Wenfei Yu, Frederick K Lamb, Rob Fender, and Michiel 
van der Klis. Peak Luminosities of the Hard States 
of GX 3394: Implications for the Accretion Geometry, 
Disk Mass, and Black Hole Mass. The Astrophysical 
Journal, 663(2):1309-1314, July 2007. 

[160] Wenfei Yu and Joshua Dolence. A HardtoSoft State 
Transition during a Luminosity Decline of Aquila XI. 
The Astrophysical Journal, 667(2): 1043-1047, October 
2007. 

[161] Wenfei Yu and Zhen Yan. State Transitions in Bright 

Galactic X-Ray Binaries: Luminosities Span By Two 
Orders of Magnitude. The Astrophysical Journal, 
701 (2): 1940-1957, August 2009. 

[162] Jing Tang, Wen-Fei Yu, and Zhen Yan. RXTE /ASM 
and Swift /BAT observations of spectral transitions in 
bright X-ray binaries in 20052010. Research in Astron- 
omy and Astrophysics, ll(4):434-444, April 2011. 

[163] Alan P. Lightman and Douglas M. Eardley. Black Holes 
in Binary Systems: Instability of Disk Accretion. The 
Astrophysical Journal, 187:L1, January 1974. 

[164] N. I. Shakura and R. A. Sunyaev. A theory of the 
instability of disk accretion on to black holes and the 
variability of binary X-ray sources, galactic nuclei and 
quasars. Monthly Notices of the Royal Astronomical So- 
ciety, 175:613-632, 1976. 

[165] T. Piran. The role of viscosity and cooling mechanisms 
in the stability of accretion disks. The Astrophysical 
Journal, 221:652, April 1978. 

[166] Marek Gierlinski and Chris Done. Black hole accre- 
tion discs: reality confronts theory. Monthly Notices of 
the Royal Astronomical Society, 347(3) :885-894, Jan- 
uary 2004. 

[167] Andrzej A. Zdziarski and Marek Gierliski. Radiative 
Processes, Spectral States and Variability of Black-Hole 

Binaries. Progress of Theoretical Physics Supplement, 
155:99-119, 2004. 



33 



[168] Shigenobu Hirose, Julian H. Krolik, and Omer 
Blaes. RADIATION-DOMINATED DISKS ARE 
THERMALLY STABLE. The Astrophysical Journal, 
691(1):16-31, January 2009. 

[169] Sheng-Ming Zheng, Feng Yuan, Wei-Min Gu, and Ju- 
Pu Lu. Revisiting the Thermal Stability of Radiation- 
Dominated Thin Disks. The Astrophysical Journal, 
732(1):52, May 2011. 

[170] T. Bclloni, M. Mcndcz, A. R. King, M. van dcr Klis, and 
J. van Paradijs. An Unstable Central Disk in the Super- 
luminal Black Hole X-Ray Binary GRS 1915-^105. The 
Astrophysical Journal, 479(2):L145-L148, April 1997. 

[171] Sergei Nayakshin, Saul Rappaport, and Fulvio Melia. 
Timedcpendent Disk Models for the Microquasar GRS 
1915+105. The Astrophysical Journal, 535(2):798-814, 
June 2000. 

[172] Agnieszka Janiuk, Boena Czerny, and Aneta Siemigi- 
nowska. Radiation Pressure Instability Driven Variabil- 
ity in the Accreting Black Holes. The Astivphysical 
Journal 576(2):908-922, September 2002. 

[173] Da-Bin Lin, Wei-Min Gu, and Ju-Fu Lu. Effects of the 
stress evolution process on the thermal stability of thin 
accretion discs. Monthly Notices of the Royal Astronom- 
ical Society, 415(3) :2319-2322, August 2011. 

[174] Marek A Abramowicz and P. Chris Fragile. Foundations 



of Black Hole Accretion Disk Theory. Living Reviews 
m Relativity, 16(1):1, April 2013. 

[175] R S Nemmcn, M Georganopoulos, S Guiricc, E T Meyer, 
N Gehrels, and R M Sambruna. A Universal Scaling for 
the Energetics of Relativistic Jets from Black Hole Sys- 
tems. Science, 338(6113):1445-1448, December 2012. 

[176] Bing Zhang. Gamma-Ray Bursts in the Swift Era. Chi- 
nese Journal of Astronomy and Astrophysics, 7(l):l-50, 
February 2007. 

[177] Jin Zhang, En- Wei Liang, Shuang-Nan Zhang, and J. M. 
Bai. RADIATION MECHANISMS AND PHYSICAL 
PROPERTIES OF GeVTeV BL Lac OBJECTS. The 
Astrophysical Journal, 752(2): 157, June 2012. 

[178] Shuang Nan Zhang. Similar phenomena at different 
scales: black holes, the Sun, 7-ray bursts, supcrnovae, 
galaxies and gala:xy clusters. Proceedings of the Inter- 
national Astronomical Union, 2(14):14-25, December 
2007. 

[179] Yuan Liu and Shuang Nan Zhang. Dusty Torus For- 
mation By Anisotropic Radiative Pressure Feedback 
of Active Galactic Nuclei. The Astrophysical Journal, 
728(2):L44, February 2011. 

[180] Hua Feng and Roberto Soria. Ultraluminous X-ray 
sources in the Chandra and XMM-Newton era. New 
Astronomy Reviews, 55(5-6):166-183, November 2011. 



